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_ , Abstract 
(D , 

r~| , We study extremal and non-extremal generalizations of the regular non-abelian 

monopole solution of [^], interpreted in as 5-branes wrapped on a shrinking S*^. 
^ ! Naively, the low energy dynamics is pure Af = 1 supersymmetric Yang-Mills. How- 

ever, our results suggest that the scale of confinement and chiral symmetry breaking 
in the Yang-Mills theory actually coincides with the Hagedorn temperature of the little 
string theory. We find solutions with regular horizons and arbitrarily high Hawking 
temperature. Chiral symmetry is restored at high energy density, corresponding to 
large black holes. But the entropy of the black hole solutions decreases as one proceeds 
to higher temperatures, indicating that there is a thermodynamic instability and that 
the canonical ensemble is ill-defined. For certain limits of the black hole solutions, we 
exhibit explicit non-linear sigma models involving a linear dilaton. In other limits we 
find extremal non-BPS solutions which may have some relevance to string cosmology. 
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1 Introduction 



One of the main motivations of the AdS/CFT correspondence and its generahzations 0, ^ 
(see 1^ for a review) is to give an account of confinement based on string theory [^]. Since 
the duahty is most naturally formulated for strongly coupled gauge theories, this goal might 
not seem too distant; and indeed, there have been many attempts, starting with to give 
a qualitatively correct description of confinement based on semi-classical reasoning on the 
supergravity side. 

A particularly natural venue for making an explicit connection between string theory and 
gauge theory is pure J\f = 1 super- Yang-Mills model. This theory exhibits chiral symmetry 
breaking and confinement, but supersymmetry gives enough control to make a number of 
exact statements (see, e.g., for a review). In particular, for gauge group SU{N), there is 
a Z2N chiral /^-symmetry (acting as a complex phase on the gauginos) which is the remnant 
of the U{1)r of the classical theory after instanton effects are taken into account. A choice 
of vacuum breaks this further to Z2 through a gaugino condensate, (tr AA) = e^^A^, where 
k = 1, . . . ,N labels the vacua, and A is the dynamically generated scale. For high enough 
temperatures, the full Z2N chiral symmetry should be restored, and we should have (tr AA) = 
0. 

The original motivation for this paper was to study the chiral symmetry breaking tran- 
sition of A/" = 1 super- Yang-Mills theory via a supergravity dual. 

In the recent literature, there are two particularly notable attempts to provide super- 
gravity duals of A/" = 1 super- Yang-Mills theory ||^, ^.[] In the geometry is a warped 
product of R'^'^ and the deformed conifold, which is supersymmetric 0, |15|, and can be 
thought of as the result of wrapping M D5-branes on the S"^ of the conifold's base Jl^ and 
then turning on other fields |18| to keep the dilaton fixed . The S"^ shrinks, but the three- 
form Ramond-Ramond (R-R) flux from the D5-branes remains; also there is a R-R five-form 
corresponding to an indefinite number of D3-branes which grows with energy scale. The 
gauge theory dual involves a "duality cascade" of SU{N) x SU {N + M) gauge theories with 
A^ = 1 supersymmetry, where N also grows with energy scale. At low energies, only pure 
SU{M) gauge theory remains. In 21, 22| an understanding of chiral symmetry restora- 
tion at high temperature was reached: black holes were shown to exist which corresponded to 
thermal states in the gauge theory with exactly zero gaugino condensate. Unfortunately, the 
supergravity equations that determine these black holes are formidable coupled differential 



equations, and the best that could be done |]22| was to solve them in a high-temperature 
expansion. This leaves open the nature of the chiral symmetry breaking phase transition. 



""^In these papers the supergravity backgrounds have non-trivial dependence on the radial coordinate ("en- 
ergy scale") only. An earlier approach, based on a massive deformation of A/" = 4, has been studied in 
[ p^ PH , p2| |l3|. The ten-dimensional supergravity solutions are more complicated in this approach because 
there is angular as well as radial dependence. Studying finite temperature in these backgrounds is difficult; 
see however M. 
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The current paper focuses on the other approach 0, which was based on re- 
interpretation of a supergravity solution previously found in |^3|, |3 . Here the R-R five-form 
field is turned off altogether, and only the D5-branes remain. The S-dual NS5-brane version 
of this geometry (with the R-R two-form replaced by the NS-NS two-form) falls p5|, ^ into 
the general category described in [^. The gauge theory interpretation is that one starts 



with little string theory ||2^ on the six- dimensional D5-brane worldvolume and compactifies 
on 5*^ to obtain four-dimensional supersymmetric Yang- Mills theory (for a discussion of some 
properties of this theory see |l30|). 

The approaches of [|, ^ thus provide different UV completions^] of pure N = 1 super- 
Yang-Mills theory which can be studied in string theory via extensions of the AdS/CFT 
corresp ondence . 



1.1 Summary of results 

It may seem that the approach of should be simpler than the duality cascade of [Q. 
Indeed, it is technically simpler on the supergravity side, and we shall obtain results on non- 
BPS solutions which are considerably more detailed than the ones available for the duality 
cascade. However, our results suggest that the Hagedorn temperature of the little string 
theory either coincides or nearly coincides with the critical temperature for chiral symmetry 
breaking, so that the super- Yang-Mills modes are not cleanly decoupled from massive modes 
in its parent theory.^ This is a particularly sharp manifestation of a persistent problem 
observed in supergravity duals of confining gauge theories: generically there is not a clean 
separation of scales between higher-dimensional modes and gauge theory phenomena. A 
general argument that this should be so is that for supergravity to be valid, the 't Hooft 
coupling should be large, so if the extra matter fields freeze out at a scale A, then the scale of 
confinement is roughly e~'^''/*^^^YM)A, where cq is some constant of order 1. We may suspect 
that the "AdS-QCD" enterprise teaches us at least as much about the UV completions (in 
our case, little string theory on S"^) as it does about the low-energy confining gauge theories. 

Besides the intrinsic interest of little string theories, there are two reasons why the 
example of 5-branes on a two-sphere deserves further study. First, this system dual does 
exhibit chiral symmetry breaking in its supersymmetric ground state, and (as we shall see) 
possesses chiral-symmetry restored states at high energy density; so we have a reasonable 
shot at describing the interesting chiral symmetry breaking phase transition. Second, it 
is possible to quantize Dl-branes in the background under consideration, using (in S-dual 
language) nothing more than non-linear sigma model techniques. This is not quite ideal: 
"weaving together" planar graphs for the gauge bosons leads to worldsheets for funda- 
mental strings, whereas Dl-brane worldsheets are related to the dual magnetic variables, and 

^We use the term "UV completions" loosely here since M — lis already renormalizable and asymptotically 
free, so it doesn't strictly require any additional fields in the ultraviolet. 
■^We are grateful to I. Klebanov for a useful discussion of this point. 
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external magnetic charges are screened rather than confined. Still, it is a real novelty to be 
in possession of string backgrounds for a confining gauge theory which do not require R-R 
fields: from this S-dual point of view we have fundamental strings moving in S'^-wrapped 
NS5-brane background. 

In a first attempt was made to construct a non-extremal black hole generalization 
of the supersymmetric solution of [y, ^ ^ . Here we shall present a systematic study of such 
solutions, which extends the work in in several directions. Rather than working in ten 
dimensions, it is useful to go back to four by integrating over the threaded by the three- 
form fiux and also dropping the spatial R'^ factor (which is possible as long as we are only 
interested in questions about translation invariant quantities in a thermodynamic limit). The 
4-d framework allows us to be guided by intuition about structure and properties of familiar 
black-hole solutions.^ Indeed, the BPS solution arose from lifting a non-abelian gravitating 
monopole in four-dimensional A/" = 4 gauged supergravity back up to ten dimensions. This 
monopole is one of the few analytically known classical supergravity solutions involving 
both non-abelian gauge fields and gravity. For a review of such solutions, both analytic and 
numerical, see |^ . 



Our approach will be to consider black hole solutions with asymptotics similar to the 
gravitating monopole solution of For the most part our non-BPS solutions will be 



numerical. As we shall explain, unbroken chiral symmetry is equivalent to having only abelian 

gauge fields in the supergravity solution: the non-abelian gauge fields yield an order parameter 

for the transition. There is a critical value (depending on the normalization of the dilaton) 

of the entropy of a black hole solution below which non-abelian gauge fields must appear. At 

this critical value, a long throat develops in the geometry which is, in the string frame if we 

are describing NS5-branes on S'^ (or in Dl-brane frame if we are describing D5-branes), the 

two-dimensional dilaton black hole geometry times x T^'^. Here the T^'^ space has 

the same symmetries and topology as the familiar T^'^ base of the standard conifold [|^. Its 

metric is only slightly different. The throat solution at the critical value of the entropy is 

available analytically, and we are also able to provide a worldsheet sigma model description 

of it as well as a description of how it is deformed as one departs from the critical point. 

One might hope to map this "critical point" in the space of supergravity solutions to a 

second order chiral symmetry breaking transition in the Af = 1 gauge theory. This does not 

work out because the temperature of the critical point is actually higher than the Hagedorn 

temperature of the little string theory, which can be read off as the limiting temperature 

of black holes far from extremality. Rather, it seems that the black hole solutions we find 

help characterize little string theory on S*^ above its Hagedorn transition. It is not possible, 

in a near-horizon limit, to proceed to T > in classical, non-extremal, fiat NS5-solutions. 

^In practice, since we are interested in static spherically symmetric solutions, we will end up, as in |5[ 
|2^ , with 1-d effective action for the radial evolution of the unknown functions in the metric and matter 
field ansatz. 
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However, it seems that wrapping the NS5 on an S"^ changes the story and allows us to 
characterize higher temperature states without resorting to string theory corrections, as in 



35| . The specific heat is negative for the black holes we find, so that the entropy decreases 
as the temperature rises. This is reminiscent of speculations that at very high energies, 
string theory has very few degrees of freedom. The thermodynamic instability that negative 
specific implies seems likely to be refiected in tachyonic modes of the black hole solutions 
[]36| , p7| , [38| , p9| , similar to the local Gregory-Lafiamme instability. We postpone a detailed 
investigation of this point, focusing instead on translationally invariant questions. 

If all our black hole solutions describe effects in little string theory, then what, one may 
ask, describes the chiral symmetry breaking transition in field theory? There are no black hole 
solutions whose Hawking temperature is less than the Hagedorn temperature of little string 
theory. Thus, semiclassically, the solution that may be expected to dominate the path integral 
at lower temperatures is the original vacuum solution of f^, ^ , periodically identified in 
Euclidean time. This solution does have broken chiral symmetry. There are no globally 
regular solutions without horizons that have unbroken chiral symmetry. Thus the transition 
which restores chiral symmetry occurs precisely when one reaches the Hagedorn temperature 
and can form the abelian black holes. At this point it is only a question whether such 
black holes are entropically favored over the periodized vacuum solution. They in fact are, 
so we may provisionally conclude that chiral symmetry restoration and deconfinement occur 
simultaneously, at the Hagedorn temperature of the little string theory, and that the transition 
is first order.0 These results are in line with the familiar conclusion |^ that solutions with 
regular horizons describe a deconfined phase, while horizonless solutions describe a confined 
phase. We will revisit this issue in section 0: as we shall see, some refinement is necessary 
on account of the thermodynamic instability of the black hole solutions. 



1.2 Organization of the paper 

In section § we shall describe the class of ten-dimensional backgrounds we are going to 

consider. These IIB backgrounds involve only the metric, the dilaton, and a three form 

field strength, which by S-duality may be taken to be the R-R field strength or the NS 

one. The ansatz will be translationally invariant in three spatial direction as well as in the 

time direction, but generally the Lorentz group 5*0(3, 1) will be broken to 5*0(3) by non- 

extremality (that is, finite temperature). The six extra dimensions comprise a radial direction 

r and a transverse compact 5-d space with 5^ x 5^ topology and SU{2) x SU{2) isometry. The 

resulting background may be interpreted [^j as a special kind of 3-brane representing D5 (or 

^It is possible that spatially non-uniform black hole solutions may have a lower minimum Hawking tem- 
perature, in which case our conclusions would be somewhat modified. It is almost certain that spatially 
non-uniform solutions play a role in describing the high temperature phase, since the specific heat is negative 
there. 

^The conclusions of this paragraph were arrived at in discussions with I. Klebanov. 
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NS5) branes wrapped over a shrinking S'^. Our general ansatz for the supergravity fields will 
be parametrized by 9 functions of the radial coordinate r, and we will derive the effective 1-d 
action for them that reproduces the full set of supergravity equations in this case. We shall 
then consider a subset of backgrounds with only 3 independent functions which corresponds 



to the solution of |2|, |2|, |[. 

In section ^ we shall obtain the equivalent set of equations from the D = 4 perspective 



23| , p4| : by looking for non-abelian black-hole type solutions of the bosonic SU{2) x f/(l) 
sector of A/" = 4 gauged supergravity (which can be obtained by compactifying D = 10 
supergravity on x T^). We shall explain the translation between the D = 10 and D = 4 
descriptions. 

In section ^ we shall study the extremal (or "zero-temperature") solutions of this sys- 
tem - solutions which have S0{3, 1) Lorentz invariance. They are obtained when a non- 
extremality parameter a is set equal to zero. We shall first consider a subset of BPS solutions 
(section 4.1) which solve a first-order system following from a superpotential and preserve 
J\f = 1, D = 4 supersymmetry. The family of these BPS solutions is parametrized by one 
essential parameter c; solutions with generic values of < c < oo are singular non-Abelian 
backgrounds, while the boundary points of the family corresponding to c = and c = oo 
are, respectively, the regular non-Abelian and singular Abelian solutions of |2^, It is 
the regular non-Abelian solution that was interpreted in as supergravity dual of A/" = 1 
supersymmetric gauge theory. 

Non-BPS (supersymmetry-breaking) extremal solutions will be described in section 4.2. 
We shall start with two "fixed-point" abelian solutions, one of which has a remarkably simple 



world-sheet description in terms of special kind of gauged WZW model |^ and 

thus is expected to be an exact string solution to all orders in a'. We shall then describe a class 
of regular non-extremal solutions (depending on one parameter b) by analyzing asymptotics 
at r = and r = oo and interpolating between them numerically. Presumably, these solutions 
may be interpreted as "excited states" of the regular BPS solitonic background, similar to 
higher excitation modes of BPS monopoles. They may be related to supersymmetry-breaking 
deformations of A/" = 1 supersymmetric gauge theory dual to the regular BPS background. 

In section ^ we shall turn to non-extremal solutions {a ^ 0) with regular black hole 
horizons. We shall determine their short- distance behavior, which depends on the two essen- 
tial parameters {Kh, Wh), the second of which may be interpreted as the U(l) chiral symmetry 
breaking parameter. The global form of the solutions is found by numerical integration. We 
shall then compute the corresponding Hawking temperature as a function of the two hori- 
zon parameters. As we will explain, there is a minimal non-zero value of the temperature, 
Tc = l/47r, which is achieved in the limit of large black holes. For Kh < oo one has T > T^, 
and the minimal value of T for a fixed is achieved for the Abelian solution, suggesting 
restoration of chiral symmetry on the gauge theory side. The limit T — » oo will lead to 
globally regular non-Abelian solutions, which break the chiral symmetry. 
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In section 6 we compute the energy and free energy of the black holes we have found. 
Remarkably, of the two-parameter family of black hole solutions, only a discrete series of 
one-parameter families has finite energy. Non-abelian black holes exist only with energy less 
than a certain threshold; abelian black holes exist only with energy greater than a different 
threshold — lower than the first, so that there is a range of energies where both abelian and 
non-abelian solutions are possible. 

In section 7 we shall address the question of chiral symmetry restoration at tempera- 
tures higher than the Hagedorn temperature. We compare the free energy of a black hole 
solution with the free energy of the globally regular BPS solution with the same periodicity 
in Euclidean time at infinity. The thermodynamic instability of the black holes makes it 
difficult to discuss Hawking-Page transitions meaningfully; however we describe conditions 
under which black holes would be expected to form. 

Section 8 contains a summary of different solutions we obtained and a discussion of 
possible application of excited monopole solutions in string cosmology context. 



While this paper was in preparation there appeared another discussion of a possible 
relation between non-extremal NS5 on S"^ background and issues of little string thermody- 
namics. There is some overlap with our section ^, to the extent that [49| also reached the 
conclusion that the specific heat is negative. We also make contact briefiy with the analysis 
of in section 1^. 



2 Ten-dimensional description of 5-branes on S'^ 

We shall study solutions in the following subsector of the type IIB supergravity action: 

Sio = \jd''xV^[R-\id^r-^e~^Hl-^e^F^) . (1) 

Here H3 = dB2 = \ Hmns dx^'^ A dx^ A dx^ and F3 = dC2 = \ Fmns dx^ A dx^ A dx^ . The 
line elements in the Einstein frame (used in the above action) and in the string frame are 
related by ds\Q^ = e~^^'^ds\Qg = gMNdx^'^ dx'^ . We shall be studying solutions with either 
F3 or being zero, so this is a consistent truncation of the type IIB theory.[] These two 
cases, i.e. the NS-NS and R-R backgrounds, are related by S-duality: if (gMTV, ^, H3, F3 = 0) 
is a solution of the field equations, then interchanging H3 ^ F3 and changing $ ^ — $ 
gives another solution with the same Einstein-frame metric gMN (but the string frame metric 
changes). In what follows we shall mostly consider the R-R version of the solutions. 

We shall be considering 3-brane-type solutions with 1+3 "parallel" directions (t, x") and 
6 transverse directions (r, 9i, 0i, ip, 6*2, ^2) representing a manifold with topology Hx x 
and metric similar to conifold metrics [133, UM. We shall assume that the metric and matter 



^Since the solutions we shall be discussing will have only metric, dilaton and one three-form non-trivial, 
they can be embedded into J\f = I D = 10 supergravity. 
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fields have non-trivial dependence on the radial direction r only, while all angular dependence 
will be fixed by global symmetries. 

Let (6'i,0i) be the standard coordinates on S"^, and {ip, 62,4)2) be the Euler angles on 
S^. We choose the 1-form basis on S"^ as (61,62), 

61 = d9i, 62 = — sin 9id(f)i , 63 = — cos 9id(pi , (2) 

where 63 is the spin connection, and the invariant 1-forms on as 

ei = cos^/'(i6'2 + sin?/'sin6'2(i02, 62 = — sin ?/'(i6'2 + cos?/' sin 6'2(i02, €3 = dip + cos 92d(t>2 ■ (3) 

These forms satisfy the Maurer-Cartan equation dea + \eabc ^bAec = 0. Let r be the transverse 
to the brane radial coordinate, while t and x" are the time and three longitudinal coordinates. 
We shall consider metrics of the following form 

dsloE = -Yi dt^ + ^2 rfx'rfx" + ^3 dr^ + n {e\ + el) + F5 (e? + e') + 63 , (4) 

where 

ei = ei - I7 61, £2 = 62 - Y-j 62, 63 = 63 - 63, (5) 

and Yi = Yi{r) are seven functions of the radial coordinate r only. 

Strings in such metric may describe confining gauge theories 0, provided Yi and Y2 
have finite limit for r — 0. That means one has finite fundamental string tension in the IR 
limit in dual gauge theory. 

In the "extremal" case of Yi = Y2 one has Lorentz invariance in 1+3 dimensional part, 
while non-extremal black- hole type solutions should have 1^/1^2 T^const. The regular horizon 
case should then represent finite temperature gauge theory in a deconfined state.^ 

This general class of metrics includes ||2^, as special cases all 3-brane-on-conifold 
metrics recently studied in the literature. For example, the subclass with Y^ = Y^^Yj = 
contains metrics whose transverse 6-space is the standard Ricci-fiat conifold, ds^ = dr"^ + 
r^dM^)^ where the base M5 = T^'^ = SU{2) x SU{2)/U{1) has 

dM^ = -{d9l + sin^ 9id(j)l) + -M + sin^ ^gc^^o) + + cos^irf0i + cos^2t^02)^ • (6) 
6 6 9 

Resolved conifold corresponds to 14 7^ Y5, = 0, and deformed conifold has F7 7^ 0. 

For = the metric has additional U{1) symmetry under ■?/'—>'?/' + ipQ, which should 

correspond to chiral symmetry on the gauge field theory side |^. If Yj{r) —>■ for r 00, 

this may be interpreted as a supergravity manifestation of chiral symmetry restoration in the 

high energy (UV) limit. As we shall see below, the symmetry under ip ^ ip + ipo may be 

restored also for Y7 = 1 and Yr, = Yq. 

^An alternative option for a finite temperature state is Y1/I2 =const with t replaced by periodic euclidean 
time. 
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In addition to the metric, we shall make the following ansatz for the closed R-R 3-form 

^3 m (y = f ) 



p 



es A {ei A 62 + ei A 62 - Ysi^i A 62 - 62 A Ci)} 
+ Yodr A (ei A Ci + €2 A 62 



(7) 



or, in terms of ci, £2, 



P 



is A {ei A e2 + - 2^71^8 + l)ei A 62 + (Ft - >8)(ei A 62 - £2 A d)} 

+ Y^dr A (ei A d + £2 A 62 



Here P is a constant which may be interpreted as a charge of D5-brane wrapped on S"^. Note 
that dFs = for any function Yg = Y8[r). Finally, we shall assume that the dilaton may be 
also non- const ant: $ = Yg{r). 

The global symmetries of our background allow one to derive all supergravity equations 
from a single 1-d effective action for functions Yi. Inserting the above ansatz for the metric and 
the matter fields into the action (|l|), integrating over all coordinates except r and dropping 
the surface term (and the overall volume factor) gives the effective one-dimensional action 
Si = J dr L, where 

L = Y: G.kiY)YX - U{Y) ^T-U. (9) 



i,k 



The action has the residual reparametrization invariance r — >■ r(r) unbroken by our ansatz. 
Expressing the Fj's in terms of 9 other functions x, 2;, /, w, w, $ 



Yi 



2z-6x 



Yo 



Y. 



^10y-2z+2l 



Y = g2s/-22+2p+2g y — ^'2y-2z+2p-2q 



Y7 = W. 



■5 = e 
Yr = w. 



Y. 



^2y-2z- 



Yo = 



(10) 



to make Gik diagonal, one finds (equivalent action was given in p5[ ) 



/2 



u 



ie' 



|g-i2p [g4g ^ ^-4,^^2 _ ^)2 ^ 2w\l - e'^P-^^Y] - 8e-2Pcosh 2q} 

^ p2 ^M+Ay+Ap |g4, ^ ^-Aq^^2 _ ^ ^)2 ^ ^{w - w f] 



Here /, which has no kinetic term, is a pure gauge degree of freedom reflecting remaining 
reparametrization invariance (e' plays the role of an einbein). Varying with respect to / one 
can then set it to any value as a reparametrization gauge. In the gauge 



/ = 
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the equation of motion for / takes the form of the "zero-energy" constraint T+U = 0. Another 
variable with a simple equation of motion is the function x{r): it is a "modulus" of the 1-d 
action as it does not enter the potential. In the gauge / = we get 



X 



0, 



I.e. 



X = —-a r , 
4 



a 



const > 



(12) 



The constant a is the "non-extremality" parameter (the choice of its sign is of course a 
convention): note that ^1/12 = e~^^ so that a 7^ corresponds to breaking of the S0{1,3) 
Lorentz symmetry in the parallel directions in the 10-d metric. 

As is clear from the action (1),(11), the charge P can be absorbed into a constant part 
of the dilaton, and so we shall assume below that P = 1. 

We shall be interested in the special subclass of solutions with 



w = w 



which corresponds to the class of solutions including that of ||2^, |J, |^ . The consistency with 
the other equations then requires that 



q = 5p 



$ + 4z - 4?/ + 16j9 = 



in which case the equation of motion for z can be integrated to give 

7 = const 



11 1^ 1 

z = -y -2p + --fr , I.e. z = -$ + -7r 



where 7 is another integration constant. 

The functions in the "parallel" part of the metric are then 

13 11 
Yi = exp[- $ + (7 + -a)r] , Y2 = exp[- $ + (7 - -a)r] . 

Assuming that 7 + |a > 0, the point r = —00 is the event horizon (as we will see below, <l> 
is finite at the horizon). To have regular horizon, we must require that the scale of the fiat 
3-space factor Y2 is finite at the horizon.^ This gives the condition 



Introducing finally 



the metric becomes 



J 2 _ */2 
""*10E ~ ^ 



7 = - a . 

' 2 



s = 2y + 2p 



g = 2q = lOp 



-e^'^'dt^ + rfx"rfx" + e^'dr"^ + e^^^el + el) + el + el + 11 



(13) 



^Equivalently, after compactifying on 3 parallel directions, Y2 becomes a scalar in 7-d theory, and, in view 
of the "no-hair theorem" intuition, one would expect that 7-d black hole will have a regular horizon only if 
that scalar does not have a charge at infinity. 
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where ^ 

$ = s — fif a r , 

y 2 

while the 3-form is given by (§) with Y-j = Y^ = w = w. 

We are finally left with only three independent functions s, and w, whose dynamics 
is determined by the Lagrangian 



/2 /2 -2(7 /2 1 4s 

s q e e 



-4^;(^2_^^2_2e-2c,_^ 



f-f/. (14) 



The only effect of the integration constant a is to modify the zero-energy constraint, 

f ^IJ = -0? . (15) 
4 

3 D = 4 description: non-Abelian black holes 
in gauged A/^ = 4 supergravity 

Before we proceed to analyzing the equations of motion for the Lagrangian (p!^ , let us re- 
derive these equations using the D = 4 approach. This is motivated by the fact that the 



solution of pS, pi, PI was originally obtained in the context of the D = 4 supergravity [23 



and then was uplifted to D = 10 p4|. It turns out that the subclass of D = 10 solutions 



determined by (|13D, ( p!4|) can be obtained in a similar way - by uplifting the D = 4 solutions. 
It will be convenient in what follows to use both the D = 10 and D = 4 descriptions, and we 
shall now establish the precise correspondence between the two. 

Let us consider the bosonic part of the action of the four-dimensional A/" = 4 half- 
gauged^ SU(2)x [U(l)]^ supergravity of Freedman and Schwarz (FS) P0[] : 

- i e^^F'^^.F"'^'^ - i a *F\^F'^^'' + i e"^*) . (16) 

Apart from the gravitational field g^u, the model contains the axion a, the dilaton $, and the 
non-Abehan SU(2) gauge field A"^ with F"^^ = d^Al-d^AI^ + SabcAl^Al. The dual field tensor 
is *F^^, = ^y/—ge^i,xpF"'^'', where e^^"^^ = 1. As was shown in [^, this model can be obtained 



via dimensional reduction of the D = 10 supergravity (A/" = 1 truncation of ([I|)) on x 
(the normalizations of the kinetic terms agree after taking into account that the radius of the 
internal manifold is ^-dependent). As a result, any on-shell configuration in the FS model, 
(g^i., A^, $, a), can be uplifted to = 10 to become a solution of ten-dimensional equations 

i°The full SU(2)xSU(2) FS model contains two independent SU(2) gauge fields |§. The half-gauged model 
is obtained by setting the second field together with its coupling constant to zero. The coupling constant for 
the first gauge field in ( |l6| ) is set to \/2, while in ||2^, Q it was set to one. The full FS model can be obtained 
from the N' = 1, D = 10 supergravity by dimensional reduction on x [OS. 
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of motion for the action (|l|). The uphfted fields are obtained as follows. The D = 10 metric 
in the Einstein frame is given by 

dsloE = e*/' (e"' '^gf.udx^dx'' + B^G" + rfx"dx") , (17) 
where (a, b,c = 1, 2, 3) 

while are the invariant 1-forms on S^. The R-R 3-form is given by 

F3 = A A 6^ - 6" A F"^ - 2e^* *rfa , (18) 
Here = | F'^^^dx^ A dx'^, and the asterisk stands for the four- dimensional Hodge dual, 



*{da) = \^/—gs^vpS d^Sidx^ A dx'' A dx\ while = 0. The D = 10 dilaton is given by 
$ + ln4.|^ If the four-dimensional configuration is supersymmetric, then its D = 10 analog 
preserves the same amount of supersymmetry. 

This correspondence between D = A and D = 10 backgrounds may be useful for 
constructing solutions m. D = 10, provided one has some insight into how to solve the 
4-dimensional problem. In general, however, it is not easy to solve the equations for the 
action (TE), unless some simplifying assumptions are made. Let us assume that d/dx^ is 
the hypersurface-orthogonal Killing vector. In this case the most general 4-metric can be 
represented as 



dsl = g^udx^dx"" = -e'^+'^rft' + e^'^-^''hik{x) dx'dx'' . (19) 

We shall also assume that temporal component of the gauge field vanishes, Aq = 0. This 
implies that the field is purely magnetic, so that ^F^^^F"^'^ = 0, and one can therefore 
consistently set the axion to zero. We are now left with the 3-metric /ij^, the gauge field 
y4°, and two scalars X and $. The equations of motion for (|16|) imply that X is a harmonic 
function, 

ViV'X = 0, (20) 

where VjV* is the covariant Laplacian with respect to the 3-metric e~'^^hik- Since a harmonic 
function is necessarily unbounded, solutions with non-constant X are singular, or possibly 
have event horizons. Using (|r^),(|rBp, any on-shell configuration (hikyA^j^jX) gives rise to 
the solution in = 10: 

dsloE = e'^^^[-e^^dt^ + dx''dx'' + e-^^hikdx'dx'' + e''Q^] , 

F3 = A 02 A 6^-6" A F^ (21) 

Although this could, in principle, give new solutions in D = 10, the equations of motion for 
the general static fields {hik,A'^,(j),X) are still rather comphcated. 

Since they differ by a constant shift, and since shifting the dilaton is a symmetry, we denote both the 4d 
and lOd dilaton by the same letter (f>. 
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For this reason we now make a further simphfying assumption by demanding that the 
D = 4 system is spherically symmetric. In this case the most general 4-metric can be chosen 
in the form 



dsi 



.e2^c/t2 ^ ^-2X+2X^^2 ^ ^2g^^Q2 ^ ^^^2 q ^^2)1 ^ (22) 



where X, A, g are functions of the radial coordinate r. The components A'^ of the spheri- 
cally symmetric, purely magnetic gauge field can be read off from 

TaA"" = w (Ti d9 - T2 sin 9 d(p) - T3 cos 6 dep. (23) 

Here w = w{r) and Tq = are constant SU(2) generators {ja being Pauli matrices). The 
corresponding gauge field tensor is 

T^F" = dw^ {Tide - T2 sin ed(f)) - T^iw"^ - 1) sin OdO A d(f) . (24) 

If w{r) = then the gauge field is of the Abelian Dirac magnetic monopole type. If w{r) = ±1, 
then F" = 0, which implies that the gauge field A"- is pure gauge and, therefore, can be gauged 
away. Below we shall use the fact that the choice w = ±1 corresponds, in fact, to the vanishing 
gauge field. 

In order to derive the 4d equations of motion, it is convenient to redefine the variables 

as 

X = X + 2s + l, <l> = s - g -^X . (25) 

Since X is a harmonic function, its equation of motion is 

(X'e-')' = 0, (26) 

which gives 

X = Xo + a J e^dr , (27) 



where Xq and a are integration constants. Inserting the ansatz (^), ( psf ) into the action 
([T^), integrating and dropping the surface term, the result is S = Att J dt J dvL, where (cf. 

(0)) 



L = e"' f s'^ - - e-''^w'' - - g'A - - e^^^' 
* 2 2^/4 



-49(^2 _ ^)2 _ 2e-29 _ 1 

Varying this effective Lagrangian gives the system of radial equations 



1 

T 



+ =-a^e^. (28) 



(e-'s'Y 

2 



1 e4-+' (-e-^^(w;2 - 1)^ + 2e-^3 + 1^ ^ (29) 

{e-'-^'w'y = e^'-^'+'iw^ -l)w, (30) 

(e-V)' = e^^+'(-e-^5(?i;'-l)2 + e-29), (31) 

-4s" + 2e-^^w'^ + 2g'^ = e^^+^z (^^-'^9^^2 _ _ ^^-2g _ _ ^2^2/^ ^32) 

X' = ae^. (33) 
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The same radial equations can be obtained by inserting the ansatz (|2^) , (p3D into the general 
equations for the action (|16D. Notice that the integration constant a enters only the last two 
equations. Since the equations are invariant under / — i>/ + /o, s ^ s — Iq/2, a— ae~'", the 
actual value of a is irrelevant, what matters is whether a vanishes or not. Eq. (|32D , which is 
the "zero energy condition," is in fact the initial value constraint. It is sufficient to impose 
it on the initial (boundary) values of solutions of the independent equations (p9D-(pT|). The 
constraint generates reparameterizations r — > f(r), which is the residual gauge freedom of the 
ansatz (^2|), (p3|) . One can fix the gauge by imposing a gauge condition on the fields (s, /, w). 
For example, one can impose the gauge condition / = 0, in which case the equation for X can 
be integrated, X = Xq + ar. 

In the / = gauge the Lagrangian (|28|) coincides with the one ([T^) obtained within 
D = IQ approach. Let us also compare the uplifted fields with those given by Eqs.(§),(|13[) 
(identifying 9 = 9i, (p = (pi). Using the notation of Eq.(||) one has = wei, = we2, 
= 63, also = dw A ei, = dw A 62, F^ = (w^ — l)ei A 62- The 1-forms 0"^ are then the 
same as in (|): 

= ei = ei - wei, 6^ = 12 = 62-^62, 6^ = 63 = 63-63. 
Then the metric (^) takes the form 

dsloE = e*/2 \-e^^dt^ + dx'^dx'' + e^'+^'dr^ + e^'^dO^ + sin^ 6 dip") + hlc 



63 A 



61 A 62 + (1 - w^)ei A 62 + w'dr A (61 A 61 + 62 A 62) . (34) 



Setting again / = 0, in which case X = ar (with Xq = 0), these expressions are exactly the 
same as in (H), (p!3[). 

Summarizing, the four- dimensional solutions in the static, spherically symmetric, purely 
magnetic sector of the half-gauged FS model are equivalent to the "3-brane" backgrounds of 
Eqs.(P),(0). In what follows we shall study solutions for gravitating Yang-Mills fields in four 
dimensions described by Eqs.(^)- (|33|) , using ( ^If ) in order to construct their ten-dimensional 
3-brane analogs. 

Before starting to solve the equations of motion, let us rewrite them in another gauge, 
i.e. choice of the radial coordinate r. While the gauge / = is sometimes useful, in this gauge 
a finite vicinity of r = is mapped into an infinite region at spatial infinity, which may cause 
difficulties in numerical analysis. For that reason, we shall often use instead the gauge where 



Introducing the functions 



the metric becomes 



A = , i.e. l = -2s-X . (35) 

1/ = e^^ , R = e^ 

dsl = e'^'i-udt^ + ly-^dr^ + R'dn^). (36) 
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Introducing also another function 

Z = <^' 



the equations (p9D -(p3D take the following form in this gauge 



^„ 3«7'2 - R'2 R2 + 1 

R + + 



R z/R z/ 



u>" + (2Z + '^)w' - - 



, 


(37) 


0, 


(38) 


0, 


(39) 


0, 


(40) 


, 


(41) 


. 


(42) 



Z/R2 

2R2Z2 + 4RZR' + R'2 + R-(R' + RZ) - w'^ + ^""^ Z^]^^ - 

I 2a 

^' - Z 

The transformation (with constant d) 

r^e^\ $^<l> + rf, z/^e-^'^z/, w ^ w, R ^ R, Z ^ e^-^Z, (43) 

maps one solution {w(r), $(r), R(r), z/(r)} into another solution {^(e^'^r), ^(e^'^r) + 
d, R(e^'^r), e~'*'^z/(e^'^r)}. Note that in this gauge the constant term is absent in the con- 
straint (^OD but is present instead in the equation for u in (pp. 
Another obvious symmetry of the equations is (C=const) 

$ ^ $ + C , a ^ e^^a , (44) 

with all other functions remaining unchanged. Since a appears only in combination with 
e~2*, it can be set, when it is non-zero, to some fixed value by a constant shift of $. 

Finally, there is the symmetry with respect to translations, when argument of all func- 
tions is replaced as 

r ^ r + vq . (45) 



4 Extremal solutions 

Let us now study solutions of the above system of equations. There are two distinct cases: 



a = and a 7^ 0, where a is the integration constant in (|T^) or (|33|). In the first, "extremal," 
case 

a = , 

the D = 10 metric has S'0(l,3) Lorentz symmetry in the 3-brane directions. In the D = 4 
description one has X' = 0, so that u = e^^ = const. In view of the scaling symmetry (^3]) 
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(or simply rescaling t and r) one can assume, without loss of generality, that v = 1. Then 
the D = A metric (^) becomes 

dsl = e^'^{-de + dr^ + R^dn''). (46) 

Written in the string frame, i.e. without the e^* factor, the t — r part of the metric is thus flat. 
The resulting solutions are either globally regular (i.e. geodesically-complete) or have naked 
singularities. There is a special subset of BPS solutions preserving part of supersymmetry. 

For a 7^ the 4-d metric function v = e^^ is non-trivial, and we get black-hole type 
solutions that may have a (regular) event horizon. Such finite temperature solutions will be 
considered in the next section. 



4.1 BPS solutions 

The system of second-order equations following from (0), (|15D or (^) in the case of a = 
admits a special subset of solutions which satisfy the first-order system of equations, following 
from a superpotential W. As in many other similar cases, such BPS solutions preserve part 
of supersymmetry (see, e.g., |41|).Q 

In fact, in the present case, the corresponding first-order system was originally derived 
from the conditions for unbroken supersymmetry, i.e. for the existence of non-trivial 



m 



Killing spinors. In [25] the same system was obtained by first finding the superpotential for 



the action (|l^ . Since the existence of residual supersymmetry was already checked in |^ 



(with independent arguments given also in 0, |2^) below we shall follow this more transparent 
superpotential approach. 

Let us write the Lagrangian ( p8|) with a = in the form (^) 

L = GM^^-Uiy) , y^ = is,w,g), (47) 

where Gik = e~Miag(l, — |e~^^, — |). Direct inspection shows that the potential U can be 
represented as 

U = -G'" ^ . ^ ■ , (48) 

where the superpotential W is [p5|] 

1 



= ±-e^ye-49(w2 - 1)2 + 2e-29 (w;2 + 1) + 1 . (49) 



^^The existence of superpotential is related to a possibility to embed the effective 1-d system (||) into a 
globally-supersymmetric action. This, in turn, is related to the fact that we consider solutions of a bosonic 
system that can be embedded into locally-supersymmetric supergravity, as well as to special properties of 
the ansatz. Though highly plausible, in general, the existence of a BPS solution (i.e. a solution of 1-st order 
system) may not automatically imply that it will be preserving part of supersymmetry. 
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As a result, the Lagrangian ( ^7]) can be written as 

and this, in turn, imphes that solutions of the first order equations 

solve also the second-order system. 

Writing down the explicit form of the "Bogomol'nyi equations" ([5T|), one finds that the 
equations for g' and w' contain only g and w, and thus, taking their ratio, gives one first-order 
equation j^=f{g,w). Introducing 

u = w'^ , V = e^^ , 

this equation reads 

dv 

u{v + u-l)— + {u + l)v + {u-lf = 0. (52) 
du 



Remarkably, the substitution |23 



iu,viu)) ^ ip,ap)) : « = pV(^\ y^u) = -p^-u-l (53) 

dp 

reduces the problem to the simple Liouville equation 

^ = 2e^(^) . (54) 
dp^ 

As a result, one finds the following analytic form of the general solution of the first-order 
equations (|5TD: in the gauge (^) (i.e. / = —2s) the functions in the gauge field (p3D and in 
the 4-d metric (|46|) are 

w(r) = — ''^ + ^0 ^ ^ g2g(r) ^ 2(7- -)_ ^q) coth(r + ro + c) — w'^ir) — 1 , 

smh(r + ro + c) 

$ = s - c/ , e2[*('-)-*o] = e-5('-) sinh(r + tq + c) . (55) 

Here ro, c, and $o are the three integration constants for the three equations. Different 
choices of $o correspond to global rescalings of the solution, while tq can be absorbed by 
shifting r — s> r — Tq. 

The parameter c (which without loss of generality may be assumed to be non-negative) 
is essential, as different values of c lead to qualitatively different solutions. Setting c = tq = 
we obtain the globally regular solution, 

c = 0: w = — e^^ = 2rcothr-w2-l, e^^*"*") = e"f sinh r . (56) 
sinh r 
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Since w 7^ 0, the corresponding 4-d gauge field (|^) is non-Abelian. The r — asymptotics 
of this solution is 

^ = 1 - + 0{r^) , e^s = - - + ©(r*^) , e^^*"*") = 1 + — + O(r^) , (57) 

by y 

while the r ^ 00 asymptotics is given by eq. (^) below. Since the dilaton (string coupling) 
grows for r 00, for large r (i.e. in the UV) one is to switch from the R-R background 
(describing the IR region of the dual theory) to the S-dual NS-NS one with the same Einstein- 
frame metric ( ^l)) and the dilaton q^^'^ns+^o) = ^/re~'^. 

For c 7^ solutions have a curvature singularity at the point, where e^^ vanishes, and 
the parameter tq can be chosen so that e^^ > for r > 0.[^ For finite values of c these 
singular solutions have non-Abelian gauge field, while in the limit c ^ cxd we get w = 0, i.e. 
the gauge field becomes Abelian, 

w = , e^s = 2r , e2(*-*o) = ^ e" . (58) 

We have set tq = 1/2 and shifted $0 by an infinite constant (— c/2) to put solution into this 
form. Note that ([58| ) represents the large r asymptotics of the family of BPS solutions (|55|). 

We conclude that, as c is changed from zero to infinity, the family of BPS solutions 
connects the regular non-Abelian solution ( |56D with the Abelian solution (0). All these BPS 
solutions preserve J\f = 1, D = 4 supersymmetry. 



4.2 Non-BPS solutions 

Let us now consider other solutions of the second-order equations (|29| ) - (|33|) or (^7|) -(^) which 



do not satisfy (|^, and thus do not preserve supersymmetry. First note that the "Higgs" 
form of the potential for the gauge- field function w in (|T^ ) ,(|2^ ) implies that the equation (|3TD 
for w admits two simple "fixed-point" solutions, w = ±1, and w = 0. More general non-BPS 
solutions will not have a simple analytic form (a standard situation for non-BPS monopoles in 
gauge theories) and will be analyzed by a combination of short- and long-distance expansions 
and numerical interpolation. 



4.2.1 Vanishing gauge field {w = ±1) 

Let us set w = ±1. In the / = —2s gauge, the field equations (|2yD -(p3D reduce to 



s" + 2s'^ - 



-2g 



g" + 2s' g' = 



(59) 



-4s'2 -f- 2g'^ + 2e-^3 + 1 = 0, 



^^The existence of a 1-parameter family of BPS solutions which are singular for non-zero value of the 
parameter is similar to what happens in the case of fractional D3-branes on conifolds p6| . 



17 



As was explained above, for w = ±1 the gauge field can be gauged away. A"- = 0. As a result, 
there is no mixing between the S"^ and angles {ia = Ca) in the uplifted D = 10 background 



dsj^E = e*/' [-dt'' + rfx"dx" + dr^ + e^^iej + ej) + ej + el + el] , (60) 

<^ = s- g , Fs = e^ Ae^ Ae^ . 

The compact angular part of this is a direct product of S"^ x S^, i.e. the symmetry of this 
solution is enhanced as compared to all other solutions with w ^ ±1: it is invariant under 
SU{2) X SU{2) X SU{2). 

Using the third equation in ( pU[ ) to eliminate e~^^ from the first two, and introducing 
V = g' and u = s', the system reduces to 

v' = 2u^ -2uv-v^ u' = -v^ , (61) 

which gives 

The numerical solution of this equation will be described below. 

4.2.2 Special Abelian solution (w = 0, g = 0) and its NS-NS coset sigma model 
counterpart 

For w = 0, the gauge field ( p3D is of the Abelian Dirac magnetic monopole type. For / = —2s, 
a = the equations (pU|)-(P3p reduce to 



s" + 2s" = e-2^ - ^ e-'^ + ^ , g" + 2s' g' = e"'^ - e"^^ , (63) 

-4s'2 + 2g'^ + 2e-^^ + 1 - e'^^ = . 

This system does not seem to have a simple general solution, but there are two important 
special solutions. 

One special solution is already known - the Abelian BPS configuration ([58|). There is 
another simple but non-BPS solution representing background with g = 0, i.e. with constant 
radius of S"^. 

Indeed, g = solves the second equation in (|55D, and then the resulting solution is 
w = 0, ^ = 0, s=-^ + so, i.e. R = e^ = l, ^ = s-g = % + ^. (64) 



The 4-geometry ( ^61) is thus the direct product of and unit 5*^. This solution will be 
important in what follows, as it will play the role of an attracting fixed point for a class of 
globally regular non-BPS solutions. 

One may wonder if this non-supersymmetric solution is stable. In fact, the instability 
of the w = 0, g = solution is suggested by the "Higgs" form of the potential for w in the 1-d 
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action ([1^) , (|28|) . Indeed, using the fact that our background is static, and that the metric 
has 2-d Lorentz symmetry in the (t, r) plane, it is straightforward to generalize the equations 
([63| ) to the case of time t and r dependent perturbations near the solution (0) (note that 
linear s or linear dilaton provides a spatial friction term): 

- d'f6w + 6w" + V26w' + 6w = 0, -d^6R + 6R" + v^^R' - 2SR = 0, 6Z = -6R' . (65) 

w has "tachyonic" mass term, and thus its perturbations may grow with time, just as in the 
standard (w^ — 1)^ scalar potential case.Q Ignoring time dependence, the four basic solutions 
of (|65D are 

5w = exp( — r) , SR = exp( ^^^^ 

Because of the spatial friction term related to linear dilaton, 6w tends to zero for large r, 
oscillating infinitely many times as it decreases. 

The D = 10 form of this solution (written in S-dual form with F3 replaced by the NS-NS 
3-form H^) has very simple form: in the the string frame the background is the direct product 
of flat R^'^, radial r-direction with linear dilaton, and angular 5-space supported by 
flux. Explicitly (restoring the dependence on the 3-form charge P and changing the sign of 
the dilaton) 

dslo NS-NS = P {-dt^ + dx'^dx'' + dr^ + dM|) , dM^ = ej + ej + el + ej + el , (67) 

$NS-NS = -$ = -$o ^ ifa = ^ eg A (ei A €2 + ei A 62) . (68) 

This NS-NS background may be interpreted as a near-throat region of NS5-brane wrapped 
over the transverse S"^ in a special way that breaks all supersymmetries. As in other NS5 
brane cases (like the regular BPS solution (0), this NS-NS description is valid for r ^ 
when the coupling is small, while for small r one needs to consider the S-dual background 



43|. 



Like the throat region of the standard NS5-brane [0 described by R^^^ x 5^ or SU{2) 
WZW model with linear dilaton, this model has a remarkably simple world-sheet conformal 
sigma model interpretation. 

Indeed, the metric 

dMl = del + sin^ eid(j)l + del + sin^ ^2^02 + (# + cos eid(j)i + cos e2d(j)2f (69) 

is of the same ^'^^'^y^^^^^^'^) coset form as the T^'^ metric (^, but now the relative coefficients 
of the U{1) and S"^ factors are equal since this is not an Einstein space but rather a solution 
of the 5-d Einstein equations with the stress tensor term. We shall call this space T^'^. 



^^This argument does not contradict the expected stability of the w = Abelian BPS (supersymmetric) 
solution (^Sj): there g is non-trivial and w, g and <I> perturbations mix. 
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The 3-form 



= P{d'4j + cos 9i 
has potential {H^ = dB2) 

B2 = P[(cosei 
Combining the T^'^ metric 



!)i + cos92d(j)2) A {sm92d92 A d(f)2 - sm9id9i A 



;>! — cos 92d(f)2) Adip + cos 9i cos 6'2(i0i A 



(70) 



(71) 



with this antisymmetric 2-tensor we get the same D=5 NS- 
NS background that was discovered recently as a simplest representative in a special 



coset sigma models introduced in [H4|. As was checked in p3|, 



class of ^ = ^m2isu(2l 
the corresponding bosonic sigma model is conformally invariant in the one- and two-loop 
approximation (3- loop approximation in the world-sheet supersymmetric case), and there are 
good reasons to believe that (in a proper scheme) these backgrounds are exact NS-NS string 
solutions to all orders in a'. 



The string world-sheet action of this 



SU{2)xSU{2) 
1/(1) 



coset model is obtained as follows. 



Let {tjji,9i,(f)i) and {ip2,92,<t>2) be the Euler angles that parametrize the two SU(2) group 
manifolds. Taking the sum of the two SU(2) WZW models with equal levels k = P E Z and 
adding the current-current interaction term [H3| with the same coefficient P one finds 



in 



d^9id^9i + d^(t>id^(t>i + d^i^id^i^i + cos^i9^0i9,^i (r/- 



+ d^92d^92 + d^<P2d^<P2 + ^^^a^^'^a + cos92d^<P2du^2 iv^" + e^'' 
+ (cos^i9^0i + a^V-i) (cos ^2^.02 + 5.^2) {V^" + e^'^ 



(72) 



The U(l) gauge invariance of this action allows one to set '?/;2 = as a gauge choice. Denoting 
then ip = ipi, the coset model (|7^ ) becomes the same as the string sigma model corresponding 
to the D=5 target space (^) , ([H]) . 

The exact central charge of (world-sheet supersymmetric version of) this model is 



2 X 



3k' 



5 - 
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k' = k-2 , k = P. 



(73) 



k' + 2 k ' 

As in the case of the NS5 throat model, the central charge deficit of this coset model is canceled 
by the linear dilaton in (|68|). Indeed, the central charge (dilaton /3-function) equation 



lp_10) + a' 



+ 0{a 



/4\ 



(74) 

vanishes for the D = 10 background (§7^ (here D = 10, {d^^ = i = j^Hl = ^R). 

It is possible to check directly (e.g., following the discussion in [Q) that this so- 
lution breaks all supersymmetries (all such coset models in [^] were claimed to be non- 
sup ersymmetric). It may have a relation to some non-supersymmetric deformation of D=6 
little string model compactified on S'^. Returning back to the S-dual R-R background sup- 
ported by the 3-form F3, one may write down the corresponding string-frame metric as 
{9s 



"■'10 R-R 



-dt^ + rfx^rfx'^ + dr'^ + dMl) 



dp^ + -p'^i-df + dx'^dx'' + dMl) . (75) 
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One may speculate that string theory in this simple background may be dual to a non- 
super symmetric deformation of AT = 1 supersymmetric theory discussed in .|^ 

This NS-NS (or R-R) solution admits a trivial non-extremal generalization (to be dis- 
cussed below): one is simply to replace the {t,r) part of the metric and the dilaton by the 
2-d dilatonic black hole background E5|. 



4.2.3 Globally regular solutions 

Consider now general extremal non-BPS solutions of the second order field equations with 
non-constant w. For a = (i.e. u = const = 1) the independent field equations (|57|)-( 
reduce to 

R" + ^"^''^ ^ ^ - + R - 6ZR' = 0, 

^, R"-2w'2-l ZR' . 
Z' + AZ^ H h 6 



R2 R 



iw"^ - l)w 



w" + 2Zw = ' Z = <^' , (76) 

R^ 



plus the constraint 



2R2z2 + 4RZR' + R'2 - 1 - w'^ + ^ ^ - -R^ = 0. (77) 

2R 2 

We will be interested in solutions that are globally regular. This means that either the 
curvature is everywhere bounded or it takes an infinite geodesic time to reach the region with 
unbounded curvature - the spacetime manifold is geodesically complete. First of all, we shall 
consider solutions that have a regular origin, which is the point r = rg where R vanishes but 
the curvature is bounded. One can set tq = 0. The manifold cannot be analytically continued 
towards negative r in this case, and so one can assume without loss of generality that r > 0.|^ 
The inspection of the field equations shows that such solutions form a one-parameter family, 
with the following small r Taylor expansion: 



w = l-br^ + O(r^), Z = $' = 2(6^ + -L)r + 0{r^), 

R = = r - (6^ + —y + O(r^), $ = $(0) + {b"^ + —y + O(r^) 
36 12 

Here b and $(0) are free parameters. The value 



(7J 



^^While the string coupling c* = Qs^'^^ decreases for small r, as in the near-horizon D5 brane case the 
curvature grows indefinitely at r — > — oo and thus the supergravity approximation breaks down there. There 
is also the usual problem of non-decoupling (at supergravity level) of KK modes corresponding to M5 space 
since its scale is naturally of the order of the string scale. 

^^Not all globally regular solutions considered below will have a regular origin, and so the restriction r > 
will not always apply. 
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corresponds to the regular BPS solution (|56|) , while for 6 7^ | we obtain its regular, non-BPS 
deformations. Expansions (^) determine only local solutions for small r, and the next step 
is to extend these solutions to finite values of r. Our strategy will be to numerically integrate 
Eqs. ([76|) in the interval r G [0, 00) using (|78D as the boundary conditions at r = 0. Since 
the constraint (^) is fulfilled by the initial values ([78|) , it holds for all r. 

Let us discuss the boundary conditions at r = 00. Having in mind future applications, 
let us consider the general equations (|37|)-(^2D with a 7^ 0. Assuming that R — 00 for large 
r, we find the following series solutions in the vicinity of r = oo:|^ 

1 ST^ — in 9 

R = V2^-^^(l---- - + ...) + V2Vxe-^{l + - + ...) + 0{e-'n, 

V2x3/2 Ax X 

r 1 ^T^ 9 9T2 — 7 1 

$ = $oo + ^-^lna: + ^(l-^-^i - + ...)- + - + ...) + 0{e-'n , 

2 4 loa;^ ox x 

w = ^(l + ^^^ + ...)+Ca;e-"(l + ...) + 0(e-2"), 

^yx X 

- = ^-^e-^-2*-(l + ...) + 0(e-2^), x = /i(r + r„o) . (79) 



Here /i, Too, V, $00, T", C are six integration constants. Notice that 6 is the maximal number 
of integration constants a solution can have: Eqs.(|37D-(i2) can be reformulated as a system 
of 7 first order equations plus one constraint. As a result, ( [79D determines asymptotics of 
a generic solution for which R — 00 for r ^ 00. There are also solutions for which R is 
bounded for large r. It is worth noting that solutions with asymptotics ([79| ) are geodesically 
complete for large r, and moreover all curvature invariants determined by (^) vanish for 
r ^ 00. 

The parameter /i (which may be interpreted as the dilaton charge at infinity) reflects 
the scaling symmetry (^) of the equations. Comparing with (|56D , (|58D , we conclude that 
for large r the solutions generically have the same asymptotics as the BPS solutions (|5BD, up 
to a rescaling and a shift, plus the polynomial terms proportional to T, and plus also the 
exponentially small terms proportional to V. 

In the extremal case of a = the solutions for r — > 00 are then given by ( |75D with 
= 1 (we are assuming u = uq = 1). The next step is to numerically interpolate between the 



r — » asymptotics ([78|) and these large r asymptotics, to find the one-parameter family of 

regular solutions in the whole interval [0, 00). It turns out that for any 6 G (0, 1/2) the local 

regular solution (ffSD can be extended all the way up to the infinity to meet the asymptotic 

solution ( [79D with certain special 6-dependent values of the parameters /i(6), Tooip), V{h), 

$oo(&), T(6),C(6) (see Figs.|g,|g). 

The behavior of the solutions is illustrated in Fig.|l| and Fig.^ For < 6 < 1/6 the 

function w is always positive, while for 6 > 1/6 it has at least one zero. As h tends to 1/2, w 

^^These expressions apply only to solutions for which R is unbounded. Similar expansions exist for solutions 
where R is bounded. 
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Figure 1: Non-BPS solution for h = 0.2. 



Figure 2: Non-BPS solution for h 0.499. 



develops more and more oscillations around zero, while the functions R and Z start oscillating 
around their constant values (1 and respectively) corresponding to the special Abelian 

solution (|6^) . Thus, one may say that the solution (|6^) acts as the large-r attracting fixed 
point for these regular solutions. Specifically, among the four independent linear fiuctuation 
modes (|66|) near this special solution there are three modes that are regular for large r. 
These modes parameterize the "stable manifold" in the vicinity of the fixed point, and their 
existence is the reason why the nearby phase trajectories approach the fixed point. As a 
result, the trajectory that starts from the origin gets attracted by the fixed point (^) and 
stays longer and longer in its vicinity as h tends to 1/2. However, for 6 < 1/2, the trajectory 
finally gets repelled from the fixed point due to the existence of the fourth, unstable, mode 
in (|66|), and after that it goes to the region where R is infinite. 



4.2.4 Limiting solutions 

A very interesting phenomenon occurs for the special case of 6 = 1/2. For 6 — 1/2 the 
trajectory approaches the fixed point ( |64|) closer and closer, and finally for h = 1/2 the 
limiting trajectory splits into two parts. For the first, interior part the trajectory starts from 
the origin at r = 0, and in the limit r oo arrives exactly at the fixed point (|6^) - after 
infinitely many oscillations. The second, exterior part of the limiting trajectory corresponds 
to the solution that interpolated between the fixed point (0) and infinity. 

Let us construct first the interior limiting solution. Returning back to the Lagrangian 
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Figure 3: Interior limiting solution 



Figure 4: Exterior limiting solution 



p8|), we introduce the new variables p{r), f{r) related to w{r), g{r) via 



w 



coshp cos/, 
The Lagrangian then becomes 



cosh p sin / . 



L = e-M 



f 



/2 



tanh p p 



/2' 



As+l 



( tanh^p 



a 



(80) 



(81) 



2 sin^ / 2 sin^ / / 4 V sin^ / sin^ f ) 4 ' 

The advantage of such a parameterization is that, as one can immediately see, p{r) = is a 
solution of the equations of motion. This means that the field equations admit the following 
first integral 

w;2 + e^^^ = 1 . (82) 

It turns out that for 6 = 1/2 this condition arises automatically. Indeed, the equation for 
p{r) derived from (^Tj) shows that for 6 = 1/2 the function p and all its derivatives at r = 
vanish. As a result, we have p = 0, and the Lagrangian (|8TD becomes simply 



e-'f2 e- 
1 



e-'s'2 



2 sin^ / 2 sin^ / 
The field equations are then (in the gauge / = —2s) 
1 



s +2s 
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f" + 2s'f = (1 + n cot / , 2 sin' f = f" + 1 . 



^3) 



4) 



• 2 f ' 

sm / 

The solution will be regular at the origin if s = Inr + O(r^) and f = r + O(r^) for r — > 0. 
Integrating (|8^ with these boundary conditions shows that / — > 7r/2 for large r. Recon- 
structing ly, R, and Z = s' — f cot /, finally gives the solution shown in Fig.^ This solution 
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is globally regular (regular at r = 0) and for large r it tends to the special Abelian solution 

Consider now the exterior limiting solution. Here R never vanishes, so the range of r is 
to be taken from —oo to +00. The solution starts from the special Abelian solution (|6^) at 
r = —00. Eg . (1661) shows that there is only one mode around this solution which is stable for 



r = —00: 5w = 0, 5R = exp(— i^^r), 5Z = —5R'. This shows that we must keep w = 
for all r, while R, Z can deviate from the values determined by the solution ([6^), so that for 
r —00 

R = l + exp(^^(r-ro)) + ... , Z = + 1^ exp(^^ (r -ro)) + . . . . (85) 

Here tq is an arbitrary parameter corresponding to the possibility of global translations. 
Integrating the field equations with such boundary conditions shows that for r +00 the 
solution follows the asymptotic behavior (|79|); see Fig.^ 

To recapitulate, both the interior and exterior limiting solutions shown in Fig.|^ and 
Fig.^ are globally regular. The interior solution interpolates in the interval [0, 00) between 
the regular origin and the special Abelian solution (|64D. The exterior solution interpolates 
for r G {—CO, +00) between the solution ( |64|) and the asymptotic (BPS) solution ([79|). 

Summarizing this section, globally regular solutions exist for b G [0, 1/2]. The solution 
with 6 = has not been described so far: in this case w{r) = 1, which corresponds to the 
case described by Eq. (|6^) . The qualitative behavior of R and Z is then similar to that shown 
in Fig.|I]. If 6 < then solutions are still regular at the origin, but w diverges at some finite 
r, where these solutions develop a curvature singularity. For b> 1/2 solutions have compact 
spatial sections, since R develops a second zero (in addition to the one at r = 0) at some 
finite r, where the geometry is singular. This type of behavior is somewhat similar to what 
is shown in Fig.^ for black holes. 

As we shall see below, among all globally regular solutions described above, there is 
only a discrete subset of solutions for which the energy is finite. 



5 Non-Extremal solutions: Black holes 
5.1 Solutions with regular horizon 

We shall now turn to non-extremal solutions that have a non-constant function X in the 10-d 
metric (|3^ ) or u in the 4-d metric (|36|), corresponding to the case of non-zero non-extremality 
parameter a in (|33|) or (41). Such solutions generalize the regular extreme solutions described 



in the previous section to the case when an event horizon is present. Since a enters ( ^if ) in 
combination 2ae~^*, it can be rescaled by shifting $ by a constant. In particular, one can 
set a = 1/2, which we shall assume in our numerical analysis. Since u = e^^ is non-constant. 
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such non-extremal solutions may have a regular event horizon. A solution has a regular event 
horizon if there is a point r = where v has a simple zero, while all other functions are finite 
and differentiable at this point. 

Without loss of generality one can set r/^ = (since the equations are autonomous). 
The field equations then admit, in the vicinity of r = 0, local solutions characterized by the 
following Taylor expansions: 

2ae-2*'' 



r + 0{r^) , w = Wh + {2ae-^^^)-^Wh{wl - 1) r + 0{r^) 



R = R, + (2ae'^^")'' ^' r + 0{r') 

R/i 



2RI 



$ = ^h + i2ae-"^'^)-' ^ ^r + 0(r^). (86) 



The parameter k, = 2ae~^*'' may be interpreted as a characteristic "mass scale" of black hole. 
The free parameters R^, and Wh determine the value of the dilaton at the horizon, the 
"radius" of the horizon, and the value of w at the horizon. One may check that all curvature 
invariants are finite at the horizon. 

We now numerically integrate Eqs.(P7|)-(^2D towards large r using ( |5B| ) as initial values 
at r = 0. For each set of values of Wh, and Rh this gives us a black hole solution living 
in the interval r G [0,r^,], where can be either finite or infinite. The set of black hole 
solutions is therefore three dimensional and has one dimension more as compared to the 
regular solutions described in the previous section, where we had only two parameters - b 
and $(0) in (|78|). The additional parameter arising in the black hole case determines the 
radius of the even horizon. 

In order to qualitatively describe these black hole solutions for different values of Wh, 
and R/i, we first notice that choosing different values of $/i leads merely to global rescalings 
of the configurations. For this reason we can set $ft = 0, since for other values of the 
structure of solutions is qualitatively similar. 

Since the equations (|37D-(^2D are symmetric under w —w, one can assume that 
Wh > 0, and then one can show that Wh must belong to the interval [0, 1], since otherwise w 
diverges at some finite r. 

Setting Wh = 0, we will obtain Abelian solutions with w = 0, while Wh 7^ will give 
non-Abelian solutions. They are qualitatively similar, the only difference is that for Abelian 
solutions w = everywhere, while for non-Abelian ones w starts from a finite value at the 
horizon and then approaches zero for large r. As was discussed above, configurations with 
w = respect the U{1) symmetry {ip ip+ipo), so Wh may be regarded as an order parameter 
for chiral symmetry breaking. 



The horizon value of R - the parameter Kh plays a crucial role. For Kh > y 1 ~ 
the solution has the asymptotic form ([79|), such that R — > cxo for r — > 00. A typical solution 
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Figure 5: Black holes with R/j > ^1 - w^. Figure 6: Black holes with < ^1 - w^. 
This corresponds to Figure Ha. This corresponds to Figure 0b. 



of this form is illustrated in Fig.^. For R/j < y 1 — wl, the event horizon is still regular, but 
the asymptotics change completely. R is no longer unbounded, but reaches a maximal value 
at some finite r; after that it decreases and finally vanishes at some r = r^, where there is a 
curvature singularity. Such a solution is illustrated in Fig.F 



In the "intermediate" case, i.e. for R/^ = y 1 — wf^^ the function R tends, for large r, 
to a constant Rqo- The whole configuration asymptotically approaches the (rescaled) special 
Abelian solution so that w oscillates, w ~ e~^°°'' sin(Zoo(r — tq)), and Z — Z^o ~ 



R — Roo ~ e Such a solution is illustrated in Fig. 0. For Wh = Q and R/j = 1 the solution 
is easy to find analytically by solving (p^)-(|42|): 



w = 0, R = 1, Z = const, $ = $0 + ^r, p = ^-^ ^-2<^,-2Zr _ ^g^^ 

For a = 0, choosing Z = we get the extremal solution (p^). In the case of a 7^ the 4-d 
metric (^) is simply the direct product of S'^ and the 2-d dilatonic black hole background 
(with the "cigar" metric in euclidean signature case) 



For 7^ and R/^ = y 1 — the non-abelian component of the gauge field is turned 
on, leading to more general solutions which may be thought of as finite deformations of the 
"cigar" . 

The results of the previous paragraph were discovered numerically, although it may be 
possible to prove them directly by qualitative analysis of the system of differential equations. 
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Solutions of this type are sometimes called "bags of gold.' 
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Figure 7: A typical solution with Rh = Figure 8: Non-extremality a for black hole so- 



y^l — w^. The oscillations in w are matched by lutions with Rh > ^1 — normalized such that 

oscillations in R, too small to be seen in this fig- $/i = 0, ^{oo) = 1. The region above (below) the 

ure. These oscillations are depicted in Figure ^ curves corresponds to values of a for solutions 

in magnified form. with $/i > {^h < 0). 



To support the claim that for R/j = y 1 — the solution is asymptotic to the cigar geometry 
for large r, recall the parametrization (^). Putting Rh = \/^~'^h amounts to setting 



the function p in ( pO|) to zero at the horizon, and, as we saw before, this implies p = 
everywhere, so that + = 1. Linearizing the analytic solution (pT]) around w = 0, one 
finds the claimed damped oscillatory behavior, which is actually the same as in Eq.(^), so 
this solution is a stable attractor as one proceeds to large r. It turns out (as is confirmed by 



numerical analysis) that for all Wh in the interval (0, 1), R/i = y 1 — wl leads to this attractor 
at large r. A summary of the resulting picture is given in Figures ^ ^ and |^. 



One may regard the behavior as one crosses from R^, > \Jl — wl to Rh < \Jl — w\ as 
some kind of phase transition, with R^ being the order parameter. 

Having qualitatively characterized black holes in the theory, we would like now to choose 
a suitable normalization for solutions whose asymptotic behavior for large r is given by ([791). 
So far we have assumed that a = \ and $/i = 0; this choice leads to an asymptotic value 
of the metric function v which is not generically equal to one, z/(oo) ^ 1. We now wish to 
rescale all solutions in such a way that 

uioo) = 1. (88) 



At the same time, we would like to keep the value of the dilaton at the horizon fixed, since it 
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f) and g) Globally regular 
extremal solutions 

Figure 9: A qualitatively correct depiction of the "phase diagram" of black hole solutions, and of particular 
solutions. Quantitatively correct plots showing some of the same information can be found in Figsj|, 0, ||, 
p!6| , and ^ Left: In the unshaded region, solutions are asymptotic to (|79|); in the shaded region, solutions are 
singular at finite r; and on the semi-circular border between I and II, solutions are asymptotic to the cigar 
geometry. The dark lines represent those solutions for which T = in (^9|) , which means that the asymptotics 
at infinity is asymptotically close to the BPS solution. Right: The {t,r) parts of the metrics, in Euclidean 
signature, are the surfaces of revolution of the curves shown. 

determines the coupling constant on the gauge theory side. Let us assume again that = 0. 
In order to be able to fulfill these two conditions at the same time, it is necessary to allow for 
arbitrary values of the non-extremality parameter a. The procedure is then as follows. Given 



a solution with $/j = and a = | for some Wh and R/^ > y 1 — w^, for which u asymptotically 
approaches some value z/(cxd), we apply the scale transformation (|43D with d = ;^ ln(z/(cxD)). 



This maps the solution to another black hole solution for which u asymptotically tends to 
one. For this new solution we still have a = |, but $/i is not longer zero but rather = d. In 
order to restore the original value of $/i we apply the scale transformation (0) with C = —d. 
This preserves the asymptotic value of z/, but changes the value of a = | to a = |e~^'^. As 
a result, the non-extremality parameter a is now fine-tuned in such a way that we have a 
black hole solution with both $/i = and z/(oo) = 1. In Fig.|] we show the values of the 
non-extremality in such normalization for both abelian and non-abelian black hole solutions. 

In order to obtain solutions with 1/(00) = 1 and for some other value of dilaton at the 
horizon, we apply the scale transformation (|^) with C = which multiplies the vertical 
coordinate of the curves in Fig.RI by e^*'' . It follows then that for solution with > the 
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values of a belong to the region above the curves in Fig.|, while for those with $ft < 0, a is 
in the region below the curves. 



5.2 Hawking temperature 

Let us compute the Hawking temperature. Switching to the NS-NS description and passing 
to the string frame, the 10-d metric becomes (cf. 



dsl^s = -^df + dx^'dx'' + u-^dr^ + e'^^ide^ + sin^ 9 dcj)"^) + i^l^ . (89) 
Let us examine the (t, r) part of the metric analytically continued to the Euclidean region: 

dsl = u{r)dT^ + u~\r)dr^ . (90) 
Near r = we have u ~ z/V, where u' can be read off from (|86D : u' = 2Q;e~^*''/R^. As 



a result, ds"^ = v'rdT"^ + Introducing p = yAr/u' and = |z^V, the metric becomes 
ds"^ = p^dd'^ + dp^. Since d should be periodic with the period 27r, r should be periodic with 
the period 13 = Att/u', which determines the inverse temperature. In the normalization (|88D 
the metric (^) is asymptotically flat, and evaluating the temperature at infinity then gives 
= f3. If one uses some other normalization of solutions, then the temperature at infinity 
will include the additional correction factor 1 / ^Juipo), which finally gives 

T = t: , . 91 

It is worth noting that this expression is invariant with respect to the scale transformations 
(|I5|), and so it does not, in fact, depend on value of i/(oo). In addition, the temperature 
is invariant also under (|4^), and this implies that it does not depend on as well. As a 
result, the temperature depends only on the two essential parameters: T = T{wh,^h)- Here 
Wh and Rh must belong to the physical region, —l<Wh<l, \Jl — w\ < R^; this is the 
unshaded region in Fig.|^. For = 0, R/^ = 1 we have the exact solution (^), for which 
T(0, 1) = v^/47r. The numerical evaluation reveals that for a fixed R/i > 1 the function 
T{wh, Rfe) reaches its minimum for Wh = and maximum for w/j = 1. For R/^ — ^ oo the 



temperature tends to a constant value, while for R/j — > the temperature diverges; see Fig.[TO 
and Fig.|ni 

The limit R^ corresponds to the lower corners of the unshaded region in Fig.^ and 
so it requires that Wh ~^ ±1- Solutions obtained in this limit can be viewed as the globally 
regular extremal configurations of section 4, but containing in addition a small black hole in 
the center. In the limit R/i — > the size of this black hole shrinks to zero, and outside the 
event horizon the configuration tends to the globally regular solution. 



Such a phenomenon is actually well known in the theory of hairy black holes : gravi- 



tating solitons are often capable of containing a small black hole inside. The regular solutions 
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Figure 10: 4:'7iT{wh,^h) fixed R/i > 1 and Figure 11: 4:nT{wh,Rh) for abelian {w = 
for + w\ = 1. 0) solutions. 

in our case befong to a family labeled by 6 G [0, 1/2] (with BPS solution corresponding to 
= |), and which member of this family emerges in the limit ^ depends on how the 
limit is taken. For example, if we take the limit along the left or right boundary of the un- 
shaded region in Fig.|^, that is keeping Wh = ±1, then the result will be the regular solution 
with 6 = 0, i.e. with zero gauge field. If we take the limit along the circle + w\ = 1, 
then the result will the limiting solution with b = 1/2. All other possibilities lead to regular 
solutions with < 6 < 1/2. 

It is important to emphasize that the black hole configurations tend to the regular 
ones for R/i — pointwise but not uniformly, and the limit is actually singular - since it is 
eventually accompanied by the topology change. As a result, the temperature diverges in the 
limit. This is very similar to the situation with the ordinary Schwarzschild black hole with 
vanishing mass, M ^ 0, in which case the metric tends pointwise to the fiat metric, but the 
temperature T ^ 1/M —>■ oo. 

Summarizing, for all solutions in the lowest corners of the unshaded region in Fig.|^ the 
temperature diverges. In particular, one can show that if the parameters belong to the circle 
Rl + wl = 1, then 

An hm - wl T{w,, ^l-wl) = l. (92) 

Let us consider now the opposite limit of large black holes, having R/j — oo. For asymp- 
totically fiat black holes the temperature would vanish in this limit. This does not happen 
in our case since large black holes are sensitive to the asymptotic structure of spacetime. 
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while metrics under consideration are not asymptotically flat. In turns out that T{wh, R/i) 
decreases for large R^, but does not vanish and tends to a finite limit independent of Wh- 

lim T{wh, ^h) = ^. (93) 
Rh^oo 47r 

This is a numerical result, but one can show directly that the limit exists. For large 
the function R > R/j is also large, and we can expand equations (p7|)-(^2D, keeping only the 



leading terms in R. The gauge field then decouples, while the resulting equations become 

R" - ^ + - - -(R' + 2RZ) - ARZ^ - 6ZR' = 0, (94) 
R u v 

Z' + AZ' + ^-- + ^iR' + 2ZR) + 6^ = 0, (95) 

v' 

2R^Z^ + 4RZR' + R'2 + /?-(R' + RZ) = 0, (96) 

u 2v 

v' = 2ae-2* $' = Z . (97) 

The space of solutions of this system admits the following symmetry transformation: 

R^fcR, <l>^<l>-lnA;, Z ^ Z, u^u, (98) 

where is a constant scaling parameter. The limit R^ — > cc can then be understood as 
— > oo. Since the temperature (^TJ) is invariant under such rescalings, its limit for large Rh 
exists. In order to explain the value T = one has to solve Eqs. (|9^) - 



Summarizing: there is a minimal non-zero value of the temperature, Tc = which 
is achieved for large black holes and is the same for all solutions. For a finite radius of the 
horizon Rh < oo one has T > Tc, and there exist both Abelian and non-Abelian black holes, 
but the minimal value of T for a fixed R^ > 1 is achieved for the Abelian solution, with 
w = 0. The temperature of this Abelian solution increases from for large Rh to V2Tc 
for Rh = 1. For Rh < 1 this Abelian solution no longer exists and T > V2Tc. In the limit 
T —>■ oo solutions may again become Abelian, if the limit is taken along the boundaries of the 
unshaded region with = ±1. In this case the chiral symmetry will be restored. However, 
in most cases the limit T ^ oo will lead to globally regular non-Abelian solutions, which 
break the chiral symmetry. 



6 Free Energy 

Having obtained the extreme and non-extreme non-BPS generalizations of the BPS solutions 
described above, our goal is to consider their contribution to the thermodynamics. For this 
we need to compute the free energy. Passing to the Euclidean region, such that the 4-d metric 

m is (cf. (ioD) 

dsl = e2*(z/dr2 + u~^dr^ + R'^dQ^) , (99) 
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with the periodic time r G [0, the free energy F is defined by / = j3F. Here the Euchdean 
4-d action / (cf. (^)) consists of the volume and surface terms, 

^ '' Kd^ = I,ol+Isurf, (100) 



where (p collectively denotes all physical fields, and the volume integral is taken over a four- 
volume fl enclosed by a 3-boundary S. The surface term is determined by the extrinsic 
curvature of the boundary, K. If iV^ is the outward normal to the boundary S, then 

K = V^iV^ = -^d.iy^m). (101) 



vol 



We assume that the boundary S is defined by the condition that r is constant, whose value 
is large and is taken to infinity at the end of calculations. The unit normal to the boundary 
is A^'^ = y^e~*5^, the 3-metric induced on the boundary is dl'^ = e^*(z/(ir^ + K^dQ'^), and 

Let us consider first the volume term in the action, I^oi- As in any theory with local 
diffeomorphism invariance, the on-shell value of this term reduces to a volume integral of a 
total derivative, and so can be expressed in terms of surface integrals. Explicitly, using the 
equations of motion one obtains 

Sit Jn Sn Jn 

= ^pj c/r(z/R2e2*$')' = Jiim ^ /5(z/RV*$') . (102) 

Here the lower integration limit makes no contribution, since by assumption it corresponds 
either to the origin of the coordinate system for the regular solutions, in which case R = 0, 
or to the event horizon, = 0, for the black holes. 

Consider now the surface term in the action. One has for the extrinsic curvature 

^ = 4e"'*(v^RV*)', (103) 

which gives 

Lurf[^,n = -lp\im v^e-*(v^RV*)'. (104) 

Adding the volume and surface terms together and using the field equation R^e^*i/' = 2a, we 
finally obtain 

/[^,S] = -^/? hm z/(RV*)'- (105) 

This gives the on-shell value of the action in terms of the asymptotic values of the fields at 
infinity, the latter being described by (fTOl). 

Since for all solutions the dilaton is linearly divergent at infinity, the action turns out to 
be infinite. Therefore, we need to regularize it. For this we subtract the value of the action 
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for a reference background choosing the latter to be the regular BPS solution (^61). This 



is the natural choice, since all solutions under consideration can be viewed as excitations 
over the BPS vacuum. For the BPS solution the metric is given by ( p9D with R = Rbps, 
$ = $BPS5 and with u = 1. The asymptotic value of the temperature of the black hole 
solution should be matched properly with the temperature of the BPS solution, i.e. with the 
(inverse) periodicity of its Euclidean time. To do this in a systematic way, we shall assume 
that for both solutions the coordinate r has the same period (3, but in addition for the BPS 
solution the time is rescaled in such a way that an (a priori arbitrary) constant factor z/bps 
appears in the BPS metric, 

dsl = e^^sps (^^gpg ^ ^^2 ^ ^l^^dn^) . (106) 

In other words, Peff = P^y^BPS is the effective temperature of the BPS solution. 

We now repeat the same calculation of / as above, but since, in contrast to ( |99| ) , z^bps 
does not enter the Qrr component of the BPS metric (|106|) , the result looks slightly different. 
The volume part of the action is found to be 

I,oi[^BPS,n = ^/3v^Jlm(RV*$')BPS. (107) 

Since the unit normal to the boundary at r=const is now N'^ = e"*^^^^,^, which does not 
contain ^ z/bps, the surface term of the action is 

W^BPS, S] = -I /?v^ lim e-*«P«(R'e3*)'BPs • (108) 
Adding the two terms together and subtracting the result from the black hole action I[ip, S] 



in ( |105| ), we obtain the regularized value of the action: 

/ = /[VP, S] - /[^BPS, = /5 lim {uiRh^'^y - v^(R'e2*)'Bps} -^Pa. (109) 
The free energy is then defined[^ in a r ^ cxd limit: 

F ^ P-'I = -\ ;im {KRV*)' - v^(RV*)'Bpg} -\a. (110) 
Before the limit is taken, the matching conditions at the boundary S are to be imposed 



47| . These conditions require that the 3-geometries induced on E are the same for both 
backgrounds. Since the boundary is S = S"^ x S"^ with the induced 3-geometries dt^ = 
e^'^iudr^ + R'^dil^) and dP = e^'^'^^^ {udT"^ + Rlpgrffi^)^ respectively, these geometries will be 
the same if the following conditions 

e*R = e*BPSRBpg , e2*z/ = e2*sps^j3pg (HI) 

are satisfied on S. In addition, the values of the matter fields for the two backgrounds should 

also be matched at the boundary P7| . 

^^Alternatively, one could define first the value of the free energy at a given large r by dividing I{r) by the 
local inverse temperature fi^/iAr) and then take r oo. Since the factor y/ i>{r) approaches 1 quite fast, this 
leads to the same limiting expression for the F. 
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6.1 Energy and entropy 



According to the analysis of for stationary spacetimes admitting foliations by spacelike 
hypersurfaces (which is the case for our solutions), the regularized free energy obtained 
from the action as described above can be related to the energy via the usual thermodynamic 
equation 

F = E-ST. (112) 
Here T = 1/(3, S is the entropy, and E is the conserved ADM energy 

E = -^[ ^/\^\CK-'K,)dSr, (113) 

where the integration is over the 2-boundary S'^ of the 3-surface Ej. Here and "^Kq 
are the extrinsic curvatures of S'^ in the geometry under consideration and in the reference 
background geometry, respectively. It is assumed that both geometries induce the same 2- 
metric on S'^, and that the time coordinate is rescaled in such a way that the goo metric 
components at are also the same for both 4-geometries. In addition, it is required that 
the matter fields at the boundary agree or "agree up to a sufficiently high order" |^ . 



This definition of the ADM energy is quite general, it does not require the reference 
background to be asymptotically flat|^, and it agrees [4^ with the definition based on the 



asymptotic symmetries In particular, (|113D can be applied to our solutions, which are 



not asymptotically fiat. Let us therefore compute the energy for our solutions. We have the 
three-geometry on a hypersurface of constant time dl^ = e^*(i/~^(ir^ + K^dQ"^), while for 
the BPS solution this changes to dl^ = e^'^^^^'idr^ + Rlpgdfi^), The boundary of Sf is 
a 2-sphere of constant r in the limit where r tends to infinity. The 2-geometries induced on 
Sf° are e'^^R'^dn'^ and e'^'^^^^Rlpsd^'^ , respectively They agree if 

e*R = e*«Psj^BPS (114) 

at 5*^°°. This condition fixes the geometrical Schwarzschild radius of the boundary. The goo 
metric components for the two backgrounds agree if 

e2*iy = e2*BPSj,3pg. (115) 

In addition, the matter field functions $ and w should also agree at S^, or at least the 
mismatch should tend to zero fast enough as expands to infinity. Notice that these 
matching conditions are equivalent to those in (|111|) required in the calculation of the action. 



The unit normal to is n'' = v^e-*5^, such that = VkU^ = ^e-3*(R2e2*)', 
while for the BPS we have "^Ko = Rlpse'^^^^"^ (Rlpse^*''^'')'. Inserting this into (^) and 



^"For static 4-metrics written in Schwarzschild coordinates, ds^ = —A'^{r)dt'^ + + r'^dVi^, Eq.(113) 

reduces to i? = —Vm^r^ooTA{^/B — \/Bct), where Bq refers to the reference background. For example, for 
Schwarzschild-de Sitter solution with A = B = I - 2M/r + Ar^ and Bo = 1 + Ar^ this gives E = M. 
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taking ( p.l4| ) and (|115|) into account, gives 

E = ~ \un {uiR'e'")' - V^si^'e'")'^^,} . (116) 

This is in exact correspondence with the first term in ( |110| ), and so our calculations of the 
energy and free energy agree with each other and with the general thermodynamic relation 
( |112| ), giving the following expression for the entropy of the solutions: 

S = ^(3a = 7rRy^. (117) 



Here we have used Eq.(91) for the Hawking temperature T = 1/ j3 (assuming that z/(oo) = 1). 
Since R/je*'' is the invariant geometrical radius of the event horizon, the entropy is equal to 
a quarter of the geometrical area of the event horizon. Notice that the energy and the action 
do not change under translations of r (^), while under (0), $ ^ $ + C, a ^ e^'^a, both E 
and / acquire the overall factor e^^ . 

Let us now use the above expressions in order to evaluate the energy and free energy. Let 
us choose a non-BPS solution and shift its radial coordinate to set Too = in ([791). The BPS 
solutions actually comprise the two-parameter family. One parameter in (^) is $0) which 
represents the constant part of the dilaton. Another parameter accounts for the freedom to 
shift the origin of the radial coordinate, r — r + tq (see (^5])). These two parameters can be 
fine-tuned in order to fulfill the matching conditions. Indeed, let us fix a large but finite value 



of r, which specifies the position of the boundary. Then the condition ( |115|) can be fulfilled 
by the suitable choice of z/bps - so far this parameter has not been specified. Next, one can 
choose To and $o such that (|1 14|) is also fulfilled, and in addition 



$ = $BPS (118) 

at the boundary. As a result, we can exactly match the boundary geometries and the boundary 

value of the dilaton for the two solutions. The gauge field functions w and wbps will not, 

however, exactly match at the boundary, unless the boundary is strictly at infinity (where w 

and If BPS are equal to zero).[^ If the boundary is at finite r, there will be some boundary 

discrepancy Aw = w — wbpsj which will measure the fall-off rate with which the non-BPS 

solution approaches the BPS background. For the energy to be finite, Aw should tend to 

zero fast enough as r — oo. Otherwise the excitations over the BPS background will not be 

well-localized and their energy will be infinite. 

^^It is usually impossible to exactly match the matter fields at the boundary. For example, for a Reissner- 
Nordstrom black hole there is always a jump of the electric field £ at the boundary E, since £ ~ for 
the solution, while £ = for the reference background (flat space). However, the value of this jump tends 
to zero as the boundary recedes to infinity fast enough to ensure that fields at the boundary "agree up to a 
sufficiently high order" . Physically, this condition means that excitations over the background are sufficiently 
localized for the energy to be finite. If the boundary values of fields for the solution and for the reference 
background do not agree up to a sufficiently high order, the excitations are too spread and their energy will 
be divergent. 
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As can be seen from Eq. ([79|) , all non-BPS solutions approach the BPS asymptotic for 
large r. If the parameter T in (^) vanishes, then Aw ~ exp(— r) and the asymptotic 
values are reached exponentially fast. If T 7^ 0, then the exponential fall-off is replaced by 
polynomial fall-off. In terms of the Schwarzschild radial coordinate = Re* ~ e''/^, the 
excitations with T = behave as l/r^, while those with T 7^ decay only as inverse powers 
of Inr^. It is instructive to compare this, say, to the Schwarzschild- AdS solution, where the 
excitations decay as l/r^ and the energy is finite. One can then think that all solutions with 
T 7^ approach their asymptotics too slowly for the energy to be finite. This is confirmed by 
the direct calculation (see below): matching the boundary geometries at finite r and inserting 
the result into (|116|) gives E ~ r~^/^e^, which is divergent as r 00. 

The conclusion is that non-BPS excitations over the BPS background for which T 7^ 
are too much delocalized and have infinite energy. 



6.2 Solutions with finite energy 

Let us now study the special case of the solutions for which 

T = . 

As we shall see, the energy then turns out to be finite. Non-BPS solutions with T = 
exist, one example being the Abelian black holes with w{r) = 0. In addition, there are also 
non-Abelian solutions with T = 0. 

Let us first consider the globally regular solutions. These are parameterized by 6 G 
(0, 1/2). If 6 < 1/6, then w is everywhere positive, and therefore T > 0. For 1/6 < b < 0.48 
w has a zero for some finite r, and therefore (see ([79|) ) T < 0. As a result, there is a value 
of b in between, which is 6 = 1/6, for which T vanishes. If we continue to increase b, we find 



that for 6 > 0.48 the function w develops already two nodes (see Fig.|T2D such that T is again 
positive. This shows that T vanishes again for b ~ 0.48. The number of nodes of w increases 
as 6 — s> 1/2, which shows that there is a discrete sequence of values 6„, n = 0, 1, . . . , for which 
T(6„) = 0. One has 60 = 1/6, &i ~ 0.48, . . ., b^ = 1/2. The numerical plot for T(6) in Fig.[T| 
shows the first three zeros of this function. The remaining zeros accumulate near b = 1/2, 
where T{b) oscillates with a very small amplitude, which oscillations are too small to be seen 
in the figure. The other asymptotic parameters in (|79|) for the globally regular solutions - 



^(&), ^oo{b), $00(^)1 and (rescaled) C{b), - are shown in Fig.|T^. Notice that 'P(&) vanishes for 
6 = 1/6 and is positive for other values of b. $00(6) ^ 00 as 6 ^ 1/2. 

Summarizing, among all globally regular solutions there is an infinite discrete subset of 
solutions for which T = and the configurations approach the BPS background exponentially 
fast. These solutions describe the "well- localized" excitations over the BPS background, and 
their energy, free energy, and action turn out to be finite. The first such excitation is shown in 
Fig.|l5|. Applying the same argument, one finds also black holes with similar properties. These 
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Figure 12: w{r) for the globally reg- 
ular solutions. It has no nodes for 
b = 0.1; one zero for b = 0.3; two 
zeroes for b = 0.485, and so on. 



Figure 13: T(6) for the globally reg- 
ular solution. Zeros of this function 
at 6 = 6„ correspond to finite energy 
solutions. 



finite energy black holes exist for arbitrary values of R/^ > 0, but only for some discrete values 
of Wh- It is clear that such finite energy configuration will be giving the leading contribution 
to the path integral. 

Let us explicitly compute the energy for solutions with T = 0. Asymptotics for large r 
are obtained from (|7S]): Q 

R = V2^+V2Vre-^l + - + ...), = 1 - ^e'^'^*- + . . . , 



$ = <^^ + -r --Inr - V^/¥e-''{l + - + ...) , w = Cre'"^ + . . . (119) 
2 4 r 



where we used the global symmetries (^5D,(|45|) to set Too = and fi = 1. Asymptotics of 
the regular BPS solution ([56|) can be obtained by putting here V = a = (and C = 2) and 
re-introducing the two free parameters in (^) by arbitrary shifts of r and $ (r* = tq — |, 
= $0 + i) 



Rbps = ^J2{r + r^) + . . . , $bps = + ^(r + n) - ^ ln(r + r,) + . . . , 

wbps = (2r + 2r, + l)e-'^-'^*-5 + ... , i/bps = const . (120) 



^^As was already mentioned earlier, both the globally regular and the black hole solutions have the same 
large r asymptotics given by (|7|). The constant parameters there {^ocV, ...) are of course different in the 
two cases: in the globally regular case they depend on the two constants b and <&(0) in (78), while in the 
black hole case they depend on the three constants Kh,iUh,^h in (86). 
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Figure 14: Parameters V{b), C{b), 
Tooib), and $oo(&) for the globally reg- 
ular solutions. 



Figure 15: The BPS (6 = 1/6) solu- 
tion and its first finite energy excita- 
tion (6 = 0.4807). 



We want to evaluate the expression for E in ( p,16|) at some large but finite value of r under 
the conditions ( |114| ), ( p.l5| ), and (|118|) , which are equivalent to 

z/ = z/BPs, e*R = e*«PSRBPs, R = Rbps , (121) 
and then take the limit r oo. The first of these conditions allows us to rewrite the formula 



( 116 ) for the energy as 



Since 



one has 

E = - \im ^yu I (1 



a 



-r-2<I>o 



2v^ 

which gives upon differentiation 



+ . . .)(v^e'-+2*-)' - (^/f+7:e'^+'-*+2*.y 



(122) 



(123) 



(124) 



E = lim v^z^ ( {y/r + r,e'^+' 



1 1 1 



2 v^rTrT 



a 



+ lim ^/h' -—=e 



-r-2<J>o 



1 



re 



2v^ 



ene2*° 



The second condition in (|121|) in view of (|123|) reduces to 
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(125) 



(126) 



Using it, one can rewrite (|125|) as 



E = 2Pe2*- + + l lim ^)e''+2*-, (127) 



where we have set to zero those terms which clearly vanish in the limit. The third matching 
condition in (|121|) gives r^, = 2Pr3/2e-'- + . . .. In view of this, the last term on the right in 
(|127| ) reduces in the limit to (— Pe^*°°), such that 



E = Ve^'^°° + ^a. (128) 

This is the final result for the conserved ADM energy for non-BPS - either globally regular 
or black hole - solutions with T = 0. Since the energy is invariant under constant shifts of 
r, the same expression holds for solutions with an arbitrary Too in the asymptotics. If the 
dilaton is shifted by a constant, $ — > $ + C, then a ae^^ (see (PD), while V remains 
intact, and the energy therefore changes by the overall factor e^^ . 

The action for finite energy solutions is expressed in terms of the energy and entropy as 

I = (3E-S. (129) 

For the globally regular solution the entropy vanishes and a = 0, while f3 can be arbitrary, 
so that we get 

/regular = PPe'"'- . (130) 

For the black holes, the entropy is 5* = (3a/2, while f3 = 27ra~^R|e^*'' (see (|9ll) ), so that 

/bh = — PR^e2*'^+2*-. (131) 
a 

Under a constant shift of the dilaton, + C, P and Kh are invariant, while a — > ae^*^, 

so that the action acquires the overall factor e^*". 

Summarizing the results obtained above, the non-BPS solutions described in the pre- 
vious sections generically have infinite energy. However, among these solutions there are 
special solutions with finite energy. These form discrete sets; they have fields approaching 
their asymptotic values as exp(— r), and thus describe finite energy excitations over the BPS 
background. In terms of the geometrical Schwarzschild coordinate = Re*, the excitations 
decay is l/r^, which is why the energy is finite. 

Let us now describe these finite energy solutions in more detail. 



6.3 Globally regular solutions with finite energy 

In the globally regular case, the finite energy solutions comprise a discrete one-parameter 
family. These solutions can be conveniently labeled by the integer n = 0, 1, . . ., which is 
the number of nodes of the gauge field function w{r) (solutions with n = 0, 1 are shown in 
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Fig.|T5D. Such solutions have asymptotics (^) at the regular origin (we set $(0) = 0). At 
infinity the asymptotics are those given in (|79D with T = 0. Such boundary conditions can 
be fulfilled only for the discrete values of the parameter 6 = 6„ in (^) for which the function 
T{b) in Fig.[l^ vanishes, T(6„) = 0. The asymptotic parameters in ([79|) then also assume only 
discrete values corresponding to V{b), Too (6), $oo(^), shown in Fig.|l^ with b = bn. 

The ground state solution is the BPS one, with 6 = 1/6 and n = 0, since w does not 
oscillate. Then comes its first excitation for b = 0.4807 with n = 1, for which w has one 
zero at some finite r. Then follow higher excitations. We list the parameters of several such 
excitations in Tab.l. As one can see from this table, for all excitations the coefficient V is 
approximately the same,0 but $oo increases with n as approximately 3n. As a result, the 
energy grows rapidly, E ^ 0.2 x exp(6?7,). The limit = oo is reached for b = 1/2. As was 
discussed above, the solutions then change the topology, which costs infinite energy. 



Tab. 1. Parameters of the globally regular solutions with finite energy. 



n 


b 


V 











1/6 





(1 -ln8)/4 





-1/2 


1 


0.4807 


0.2304 


2.902 


7.7 X 10^ 


5.5.258 


2 


0.4996 


0.2295 


6.083 


4.4 X 10^ 


9.750 


3 


0.499991 


0.2294 


9.175 


2.1 X 10^ 


14.121 


oo 


0.5 


0.2294 


oo 


oo 


oo 



To summarize, the globally regular finite energy solutions are characterized by the num- 
ber = 0, 1, . . . of nodes of w. The ground state energy is zero, while for all excitations the 
energy is positive and rapidly increases with n. The action / = (3E also grows rapidly with 
n, where the inverse temperature (3 can be set to any value. As a result, for any given /3, the 
ground state solution gives the leading contribution to the path integral. The contribution 
of the excitations is highly suppressed.]^ 

6.4 Black holes with finite energy 

Let us now consider the black holes with finite energy. These are obtained by selecting from 
the set of all black holes considered in section 5 only those solutions for which T = (we 

^^Numerical values of the parameters of the solutions can be determined by the multiple shooting method. 
The accurate determination of V is, however, extremely involved, since V is the coefficient in front of the 
subleading terms which are exponentially small as compared to the other, leading terms. We used a simplified 
numerical procedure giving the value of V with ~ 20% uncertainty. As a result, the values of V and E given 
in the table are, in fact, approximate. Since it requires considerable efforts to improve these numbers, we 
postpone this for a future publication. 

^^Notice that the normalization is important. One can use ( ^ ) to rescale all solutions to set $oo = 0, and 
then the energy will be V , which is approximately the same for all excitations. 
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always assume that z/(oo) = 1). For any given value of the event horizon size R^, there are 



special values Wh{n, Rh) of the gauge field function w at the horizon, shown in Fig.|T^ which 
give rise to solutions with w ~ exp(— r) for large r; see Fig.^. For all other values of Wh 
one has w ~ l/v^ fo^' large r (Fig.|16D and the energy is infinite. The finite energy solutions 
therefore comprise a discrete series of one-parameter families: particular solutions are labeled 
by {n, R/i), where = 0, 1, 2, . . . is the number of nodes of w outside the black hole horizon, 
while Rh > 0. 
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Figure 16: Function w for n = 0, 1 fi- 
nite energy black holes with Kh = 0.5. 
For comparison, two other solutions 
are shown, for which w ~ for 
large r. 



Figure 17: Parameters Wfi{n,Rh) for 
finite energy black holes. For other 
values of Wh, w{r) tends to zero too 
slow for the energy to be finite; see 



Fig.|l6. 



For n = the set of such black holes consists of two branches. First, there are the 
abelian black holes, which exist for 1 < R/^ < cxd. 

Second, for small < Kh < 1-3 there are also non-Abelian solutions. For these w starts 
from some finite value at the horizon, and then exponentially quickly tends to zero. In the 
limit R/i — > the field configurations approach the BPS solution pointwise (in the exterior 
black hole region), and so in some sense they can be viewed as black hole generalizations of 
the BPS solution itself. As Rh increases, the value of Wh for such solutions decreases, and 
finally it vanishes for Kh = Rh(0) = 1.3, at which point the abelian and non-abelian branches 
merge. For Kh > 1.3 only the abelian solutions exist. 

There are also non-Abelian black holes with n > 0. For these w starts from some finite 
value Wh at the horizon, and then after n oscillations around zero exponentially fast tends 
to zero. The function w(r) for two such solutions with n = 0, 1 and R?, = 0.5 is shown in 
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Fig.|T6[ In the limit ^ these solutions approach pointwise the globally regular finite 
energy solutions described above. As increases, the value of Wh decreases, and finally for 
some finite R/i = R/i(^) the solutions merge with the abelian black holes, similarly to what 
happens to the n = non-abelian branch. 

Summarizing, all non-abelian solutions exist only for small values of R/i, and all of them 
merge with the abelian solution for R/^ = Ri^(n), where Ri^(n) are Rft(O) = 1.3, Rft(l) = 1.01, 
. . ., R^(oo) = 1. For R/i > 1.3 only the abehan solution exists. 



* =0 





Figure 18: Energy and action for the 



n 



0, 1 black holes. 



Figure 19: Entropy a, /?, and 
for the n = black holes. 



Having obtained the black hole solutions, we can compute their thermodynamic param- 
eters. The energy E (|128|) and the action / (|131|) for the n = 0, 1 black holes are shown in 
Fig.|TB| with the normalization $oo = for all solutions.^ For R/^ — the energy of the n-th 
non- Abelian black hole coincides with that of the n-th regular solution.[^ As Hh increases, 
the energy grows. For R/j = R/i(n) the non- Abelian solutions merge with the Abelian branch. 
The subsequent increase in R/i along the Abelian branch is accompanied by further increase 
of the energy. 

For all black hole solutions the action /(R^) is zero for R^, = 0, positive for small values 
of Rh, and negative for all large enough 

However, 



^^Since E and / depend on V , their values are determined with some uncertainty; see footnote 
the quahtative behavior of the E and / curves seems to be independent on the numerical scheme used. 

^^Notice that the energy of the regular solutions in Tab.l is given in the different normalization: $(0) = 0. 
Shifting the dilaton so that — 0, their energy will be E — V, where the values of V are given in Tab.l. 

^^The action vanishes for R;i — > because / j3E, where E is the energy of the n-th regular solutions, 
but /? = 1/T — > 0, since the black hole temperature diverges in the limit. 
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In Fig.|T9] we have shown the entropy S'(R/i), the non-extremahty a{Iih), the inverse 
temperature /?(R/i), and the value of the dilaton at the horizon $/i(R/i) for the n = black 
holes. In agreement with (p3D, one has /?(oo) = An. In addition, the behavior of the ratio 
(5E/S shown in this figure indicates that for large R/j the following equation of state holds: 

E = TS. (132) 

This agrees with the first law of thermodynamics, dE = TdS, since T = is constant for 
for large Kh- We therefore recover in the UV the standard NS5 brane thermodynamics. In 



Figs. 2C, 21 we also plot the energy and free energy against entropy for the n = black holes. 



As we can see, for large black holes F also scales linearly with S. 




Figure 20: Energy against entropy for Figure 21: Free energy versus entropy 

the n = black holes. for the n = black holes. 

The value $/i(R/j) is an important parameter, since it determines the value of the string 
coupling constant. It is therefore interesting to consider another normalization for all solu- 
tions. For example, instead of fixing the value $oo = one can fix $/i = 0. Using (0), this 
can be achieved by translating $(r) — > $(r) — ^h{Iih), where $/i(R/i) is shown in Fig.|l9|. The 
energy, action, entropy, and non-extremality a acquire then the factor exp(— 2$/i(R/i)), while 
the temperature remains invariant. 

To conclude this section, we have learned the following about the value of the action 
for globally regular and black hole solutions. The action of all globally regular solutions is 
non-negative, with the minimal (zero) value achieved for the BPS solution. For black holes, 
apart from those with small Rh, the action is negative. 
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7 Restoration of chiral symmetry for T > Tc? 

We have obtained the non-extremal generahzations of the globally regular BPS solution 
P^l , P^] , which we have reproduced as The solutions which have no singularities outside 
horizons are the original BPS solution, the globally regular non-BPS solutions, and the black 
hole solutions corresponding to the unshaded region in Fig.^ Of this two-parameter family 
of solutions, only a discrete series of one-parameter families has finite energy. These are 
the classical saddle points which make important contributions to the path integral. All the 
black hole solutions have temperature larger than the Hagedorn temperature of the little 



string theory, as illustrated in Figs. p!0|JTl| . Thus, as remarked already in section |1 . 1| , the 
solution that dominates the path integral at temperatures lower than Tc is the original BPS 
solution with periodic Euclidean time. (The contributions of the globally regular non-BPS 
solutions are exponentially suppressed since their energy density is finite and positive). The 
energy, entropy, and free energy of the periodized BPS solution are equal to zero in the 
classical supergravity approximation, which only indicates that they are less than 0{N'^). 
Loop effects, due to the broken supersymmetry in the thermal boundary conditions, would 
give rise to an 0(1) free energy. This is appropriate for the low-energy A/" = 1 gauge theory 
in its confining phase. Equally appropriate, chiral symmetry is broken in this regime. A 
deconfined phase might be expected to have restored chiral symmetry, and energy, entropy, 
and free energy of order A^^ — like our abelian black hole solutions. 

So far, the discussion is little different from that of 0, where it was argued that in 
global anti-de Sitter space, a low-temperature phase corresponding to empty AdS^ gives way 
to a high-temperature phase corresponding to Ad^s-Schwarzschild through a Hawking-Page 



transition that corresponds to deconfinement in the gauge theory. The main differences 
in the current context are 1) the putative high-temperature phase is thermodynamically 
unstable, and 2) the little string theory is believed to have an exponential growth in the 
number of states at high energy. For both of these reasons, the canonical ensemble is ill- 
defined above the Hagedorn temperature Tc, and it doesn't make sense to speak of equilibrium 
processes at controlled temperatures higher than T^. Thus, though it is tempting to identify 
the abelian black hole solutions (which do have T > T^) as a high-temperature, deconfined 
phase, with restored chiral symmetry, the truth is more complicated. 

Suppose that a system such as the one we describe (that is, NS5-branes on a shrinking 
S"^) were to come in thermal contact with a "heat bath" at a temperature Tbath > T^. Small 
black holes would form and evaporate continually. Eventually, through thermal fluctuations, 
enough energy would be concentrated in one region to make a larger black hole, with a 
temperature lower than Tbath- The subsequent evolution would suck energy continually from 
the heat bath until thermal contact ceased or the heat bath fell below Tc. In regions of high 
energy density, where Rh > 1.3, chiral symmetry would be restored because the only black 
hole solutions with high enough energy are abelian. In regions of low energy density, where 
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Rh < 1, chiral symmetry is broken because the only available black hole solutions are non- 
abelian. It is likely that the end state of the system would be spatially non-uniform along 
the NS5-brane world- volume, since the uniform state is thermodynamically unstable and this 



has been associated |36, 37] with the presence of a Gregory-Laflamme instability. 



It may be noted from Figs. |21| , p]8| that the action, I = (3F, is negative for large Rh^ but 
becomes positive for small R/j. This might be regarded as the signal for a Hawking- Page 
transition back to the periodized BPS solution at very high temperatures; however this is not 
a coherent interpretation since the canonical ensemble is still ill-defined. More physically, it 
is difficult to discuss a first order transition between two phases if one is thermodynamically 
unstable, since the unstable phase may not last long enough for the transition to take place. 

For very large entropy/energy density (corresponding to very large R/i), Buchel has 
claimed S = (3hE + a log E plus subleading corrections, with a < , which result was 
obtained assuming that the thermodynamic description applies. This is consistent with our 
result that the specific heat is negative. However, it also implies that F > 0, which is opposite 
to what we obtain in our analysis. Although we reproduce the energy-entropy relation in the 
leading order, the subleading terms are different, which probably indicates the breakdown of 
the thermodynamic description. 



8 Conclusions 

Let us enumerate the solutions we have found. In citing equation numbers, we sometimes 
refer only to asymptotics if the solutions were obtained numerically. It helps to categorize 
solutions according to whether they involve the non-abelian components of the SU{2) gauge 
field when expressed in four- dimensional terms. These components are determined in terms 
of a single function w{r), and U{1) C SU{2) is unbroken precisely if w{r) = 0. It happens 
that w vanishes for all r if it vanishes at the horizon, if there is a horizon, or if not, at the 
point where the radius of the S"^ vanishes. 

1. The regular supersymmetric solution, (^6l). This solution was found in |2|, ||]. It 
preserves four supercharges and has w ^ 0. Its ten-dimensional lift was shown in to 
represent 5-branes wrapped on a shrinking 5*^, and it was therefore conjectured that the 
supergravity geometries provided a holographic description for JV = 1, D = A super- 
Yang-Mills theory. The other solutions we obtain can be viewed as excitations of this 
regular BPS one. 



2. Singular BPS solutions, (pSl). These solutions preserve four supercharges, but they 
are unphysical because of a naked singularity where the S"^ shrinks to zero size. The 
abelian "Dirac monopole" solution, (0), is a special case of the one-parameter family, 
(^5]), which includes the BPS solution (|56|) as its only regular representative. 



46 



3. 



The vanishing gauge field solution, (|60D. This solution breaks all supersymmetry, but 
it has SU (2) x SU (2) x SU{2) global symmetry, corresponding to an internal geometry 
which is S"^ X S'^. 



The factorized abelian solution, (0). All supersymmetries are broken, but the geometry 
factorizes into a five- dimensional compact coset manifold, T^'^, and a non-compact piece 
with a linear dilaton. T^'^ has a bigger U{1) fiber than the conventional T^'^ metric, 
and the interpretation is that NS5-branes have been wrapped on the 2-cycle and then 
delocalized in the other directions. We find an explicit sigma model description of this 
geometry, valid in the weak coupling region. 

Globally regular non-BPS solutions, (|78D , (^). Superficially there is a one-parameter 
family of these solutions labeled by b, including the solution (^) as its one BPS rep- 
resentative. Of these, only a discrete series has w{r) falling off exponentially at large 
radius, which we have found to be a necessary condition for finite energy. For solutions 
very far down the discrete series, there is a long region which is nearly the factorized 
abelian solution, and it is closed off on the inside by the S"^ shrinking, and on the outside 
by asymptotics similar to the regular BPS solution. 



6. Abelian black hole solutions, (p6D with Wh = 0. These solutions exist only if the 
entropy density is large enough: they are parametrized by the horizon radius, R/^ > 1. 
For Kh = 1, we have the analytic solution (0), which is the factorized abelian solution 
cut off on the inside by a black hole horizon: that is, the standard 2-dimensional dilaton 
black hole times times T^'^. 

7. Non-abelian black hole solutions, (^) with Wh 7^ 0. Superficially there is a two- 
parameter family of solutions, including all the other solutions listed as limiting cases 
(though in some cases the relevant limit is only pointwise, not uniform in r — allowing for 
instance the asymptotics to change). However, only a discrete series of one-parameter 
families has w{r) falling off exponentially at infinity. Each of these one-parameter fam- 
ilies terminates at one end on the line of abelian solutions, and at the other end at one 
of the globally regular solutions. 

Many of the qualitative features of our results can be understood from Fig.^ Roughly 
speaking, the typical non-abelian black hole solution has some oscillations of w{r) in the 
region where it is close to the factorized abelian solution. This behavior is cut off at one end 
by the horizon and at the other by expansion of the throat into asymptotics similar to the 
BPS solution. 

The globally regular non-BPS solutions, corresponding roughly to excitations of a non- 
abelian gravitating monopole, are possibly significant to string theory cosmology. These 
solutions were constructed with 3 + 1-dimensional Poincare invar iance, but they have finite 
positive energy density as compared to the supersymmetric solution. This translates to a 
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positive contribution to the four- dimensional cosmological constant. To be more precise, sup- 
pose we had constructed a compact solution where some local region was well-approximated 
by one of our globally regular, non-BPS solutions. And suppose the moduli, like the average 
value of the dilaton, were fixed. Then the non-compact four- dimensional part of the solution 
would have to be de Sitter space, and the quantity E in Tab. 1 would translate into a cos- 
mological constant. The reason we were able to construct solution with 3 + 1-dimensional 
Poincare invariance was that the extra six dimensions were non-compact, so that gravity is 
non-dynamical. We can refine things a little further if we think in terms of a toy model where 
the effects of compactification are represented by cutting off our non-compact geometry at 
some large but finite tq. Solutions with T = in (^) have finite energy as tq oo, but 
other solutions do not. In short, we expect that upon fixing finite rc, the solutions in the 
discrete series would "broaden out" into sharp, deep valleys in a four-dimensional effective 
potential. There would be only finitely many minima, because for high excitation modes the 
nodes of w would fall outside the cutoff radius. Thus the final picture is a four- dimensional 
effective potential with many minima separated by high walls. 

So far we have assumed that moduli are stabilized, but so far in string theory this seems 
very hard to do. In the scenario of the previous paragraph, the cosmological constant would 
have very weak dependence on rc, because in the rc ^ oo limit the energy computed in 
Tab. 1 is finite. However it would depend exponentially on the dilaton, so each minimum 
would extend to a long, low valley. This is not much different from the conventional picture of 
the effective potential in heterotic string compactifications with broken supersymmetry. The 
novelty is that the supersymmetry breaking occurs as a non-BPS excitation of the internal 
geometry. 

One may imagine a cosmological scenario where, at some stage in the evolution of the 
universe, one finds local physics near the shrinking S"^ described well by an abelian black hole. 
As energy density decreases due to expansion, the system would have to find its way onto one 
of the non-abelian branches in Fig.^. Only if the system found the n = branch would it then 
relax into a supersymmetric minimum; otherwise it would "lock in" some oscillations of w{r), 
and relax to a globally regular, extremal solution with a non-zero cosmological constant and 
broken supersymmetry. Thus we have given at least a rough outline of how one might end up 
in a non-supersymmetric valley of the four-dimensional effective potential and not be able to 
tunnel into a supersymmetric solution.^ This mechanism is intrinsically non-field-theoretic 
because the Hawking temperature exceeds the Hagedorn temperature of the little string the- 
ory. We consider it plausible that the contribution to the cosmological constant would be 
small if the throat region, well-described by our non-compact solutions, were long; however 
this is a point which deserves further investigation. Various drawbacks remain, notably the 
usual question of why sparticle mass splittings are so much bigger than the cosmological 

^^One might in fact question whether a four-dimensional effective potential is a valid notion. We use it for 
lack of a better language. 
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constant. Also, one may worry that the thermodynamic instabihty will lead to unacceptably 
large spatial inhomogeneities. But it nevertheless would be fascinating to see whether the ex- 
cited monopole solutions could be embedded into a global string compactification — preferably 
one with other ingredients which fix the dilaton. 
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